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Abstract 

An  analytical  study  was  conducted  to  determine  the  Influence  of 
moisture,  temperature,  and  curvature  on  the  bifurcation  load  of 
cylindrical,  composite  panels  subjected  to  a  simple  shear  loading.  Two 
laminate  ply  orientations,  (0/45/-45/90 Jg  and  (45/-45]2g,  were  analyzed 
for  six  radii,  four  temperatures,  and  two  Initial  moisture  conditions. 
The  elght-ply  composite  panels  were  assumed  to  be  manufactured  from  a 
graphlte/epoxy ,  AS/3501-5.  To  evaluate  the  Influence  of  moisture  and 
temperature,  the  transverse  modulus,  E2»  and  shear  modulus,  Gj2» 
degraded  based  on  experimental  test  data  for  the  AS/3501-5  system.  Each 
ply  orientation,  for  a  12  In.  panel  radius,  was  evaluated  at  20 
time/ temperature  conditions  that  ranged  from  80  to  300*F,  and  moisture 
concentrations  ranging  from  a  zero  moisture  content  to  an  equilibrium 
moisture  distribution.  The  investigation  of  curvature  was  conducted 
only  for  the  [45/-45]2s  laminate  and  at  a  limited  number  of 
time/ temperature  conditions. 

The  bifurcation  loads  were  determined  using  the  STAGS-Cl  finite 
elements  shell  analysis  program.  This  analysis  used  the  prebuckled 
linear  displacement  option  to  calculate  the  bifurcation  loads.  An 
Increase  In  temperature  and  moisture  absorption  was  found  to  cause  a 
reduction  In  the  panels  bifurcation  load  ranging  from  a  maximum  of  25.6 
percent  for  the  [0/45/-45/90]g  laminate  to  34.5  percent  for  the 
t45/-45]2s  laminate  for  the  panels  with  a  12  In.  radius.  This  reduction 


In  the  bifurcation  load  Is  significantly  Influenced  by  the  change  in 
curvature  at  elevated  temperatures  and  moisture  contents  The  maxlnuD 
reduction  in  the  bifurcation  load  varied  from  34.5  percent  for  a  12  In. 

panel  radius  to  11.9  percent  for  a  10,000  in.  panel  radius  a  decrease  of 
22.6  percent  for  the  I45/-45]28  laminate. 
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THE  BUCKLING  OF  COMPOSITE 
CYLINDRICAL  PANELS  CONSIDERING 
ENVIRONMENTAL  EFFECTS 


I.  Introduction 


Background 

Composite  materials'  high  strength  to  irelght  ratio  makes  It 
particularly  well  suited  for  aircraft  structural  use.  The  desire  to 
Improve  the  efficiency  of  the  structural  design  through  the  application 
of  composites  In  order  to  tailor  the  component's  strength  and  stiffness 
to  match  the  load  and  stiffness  requirements  results  In  components  with 
unique  structural  responses.  The  USAF's  forward  swept  aircraft  design 
Is  an  excellent  example  of  the  materials'  advantage  In  that,  through 
proper  composite  "tailoring",  the  problem  of  aeroelastlc  divergence  can 
be  practically  avoided  without  the  enormous  weight  gain  associated  %rlth 
conventional  metals.  One  difficulty  In  such  applications  of  composites 
Is  the  Inability  to  use  conventional,  classical  structural  design 
practices  to  predict  the  structural  response  of  composite  laminates. 

The  fact  that  their  stiffness  can  be  oriented  to  preferred  directions 
make  them  complex  structures  to  analyze. 


In  a  semlmonocoque  design,  thin  skin  panels  along  with  fuselage 
frames  and  longerons  or  wing  skins  with  spars  and  ribs  a:  e  used  to 


achieve  an  efficient  structure.  In  such  a  design,  the  accurate 
determination  of  the  buckling  load  of  curved  skin  panels  is  necessary  to 
prevent  premature  structural  failure.  It  is  stated  In  Sef.  [1]  that, 
"Numerous  papers  concerning  the  Instability  of  laminated,  anisotropic 
plates  and  shells  can  be  found  In  the  open  literature.  Buckling  of 
curved  panels  has,  however,  received  little  attention.”  This  Is 
especially  the  case  for  shear  buckling  of  curved  panels.  Some  recent 
work  has  been  carried  out  on  curved  pane,ls  In  order  to  evaluate  their 
buckling  characterlsltcs  under  compression  loads  [2-5].  Prom  the 
standpoint  of  shear  instability  of  cylindrical  composite  panels.  Harper 
looked  at  shear  buckling  of  circular  cylindrical  shells  [6]  and  Whitney 
examined  shear  buckling  of  symmetric  and  unsymmetrlc  angle-ply, 
graphite/ epoxy  curved  plates[l]. 

Another  aspect  related  to  composite  material  research  is  the 
consideration  of  moisture  and  temperature  effects  on  the  overall 
structural  characteristics.  Snead  and  Palazotto  examined  the  effects  of 
moisture  and  temperature  on  the  Instability  of  cylindrical  composite 
panels  loaded  in  axial  compression  [7].  These  Influences,  sometimes 
referred  to  as  hygrothermal  effects,  fall  in  the  general  category  of 
environmental  influences.  Hygrothermal  effects  have  been  found  to 
significantly  degrade  the  mechanical  properties  of  most  organic  matrix 
compounds  [7-21]. 

Purpose 

The  purpose  of  this  thesis  Is  to  analytically  evaluate  the 


Instability  of  composite  cylindrical  panels  subject  to  simple  shear 
loading,  using  the  STAGS-Cl  finite  element  shell  analysis  program.  The 
influence  of  moisture,  temperature,  and  curvature  on  the  panel's 
bifurcation  load  %rlll  be  investigated  using  this  finite  element  program 

Scope 

The  cylindrical  panels  evaluated  in  this  thesis  are  12  in.  by  12 
in.,  eight-ply  laminates  made  of  graphite/epoxy.  The  panel's  thickness 
width,  and  height  are  held  constant  while  the  radius,  ply  orientation, 
laminate  temperature,  and  Initial  surface  moisture  conditions  are 
varied.  Six  radii,  tvo  ply  orientations,  and  a  simple  support  boundary 
condition  are  investigated.  During  the  investigation  of  the 
environmental  influences,  four  temperatures  and  two  surface  moisture 
conditions,  at  five  time  values,  are  evaluated. 


The  plane  stress  constitutive  relations  for  an  othotroplc  lamina 


can  be  written  as: 
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Where  the  Qij*s  are  the  reduced  stiffnesses  and  are  defined  In  terms  of 
the  engineering  constants  as: 

Qll  -  Ei/(1-»/i2*'2i) 
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where 


(3) 


In  any  other  coordinate  system  %d.th  the  fiber  axis  oriented  at  some 
angle  #  with  respect  to  the  structural  axes,  the  stresses  are: 


H 

^11  ^12  ^16 

«x 

0y 

>  - 

^12  ^22  ^26 

( 

«y 

^16  ^26  ^66 

V 

Where  the  Qij's  ate  the  transformed  reduced  stiffnesses: 

Qll  *  QnCOsV  +  2(Qi2  +  2Q65)sln2^  cosH  +  Q228ln*tf 
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To  extend  the  8tres8-8traln  relationship  to  a  multilayered  laminate, 
two  assumptions  must  be  made.  First,  the  laminate  Is  assumed  to  have 
perfect  bonds  between  the  lamina  so  that  no  Interlamlna  slipping  can 
occur.  Secondly,  the  Klrchhoff-Love  hirpothesls  applies.  This 
hypothesis  states  that  normals  to  the  mid-surface  remain  plane  and 
normal  to  that  surface  after  bending.  These  assumptions  lead  to  the 
strain-curvature  relationship  for  a  laminate. 


Where  the  o  superscript  Indicates  the  mid-surface  strains,  the  K'  s  are 
mid-surface  curvatures  and  Z  (See  Figure  2.2)  represents  a  distance  from 
the  panel’s  mid-surface. 


Figure  2.2  Geometry  of  an  N-Layered  Laminate 


For  a  cylindrical  panel  with  moderately  large  displacements  and  also 
moderately  large  rotations  of  tangents  to  its  mid-surface,  the  above 
strains  and  curvatures  are  given  by  Sanders'  kinematic  relations  without 
initial  Imperfections  [24]  as 
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Where  u,  v,  and  w  are  the  axial,  circumferential,  and  radial  components 
of  displacement,  respectively,  of  the  panel's  mid-surface.  The  <t>' 6  are 


rotational  components  and  are  expressed  in  terms  of  the  displacements 
as: 
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where  r  Is  the  panel's  radius  of  curvature. 

If  Eq.  (6)  Is  substituted  into  Eq .  (4),  the  stresses  in  the 
layer  can  be  expressed  in  terms  of  the  laminate  mid-surface  strains  and 
curvatures  as: 
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The  resultant  forces  and  moments  acting  on  a  laminate  can  now  be 
found  by  integration  of  the  above  stresses  in  each  lamina  through  the 
laminate  thickness.  Denoting  the  in-plane  ply  stress  by{7^>  then 


t/2 
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where  and  are  the  force  and  moment  per  unit  length  (%ri.dth)  of  the 
cross-section  of  the  laminite  as  shown  in  Figure  2.3. 


Figure  2.3  Forces  and  Moments  on  a  Laminate 


The  forces  and  moments  sho%m  in  Figure  2.3  can  be  written  as: 


where  N  is  the  total  number  of  laminae,  and  and  Zj^_j  are  defined  in 
Figure  2.2. 

The  integration  Indicated  by  Eq's.  (11)  and  (12)  can  be  rearranged 
to  take  advantage  of  the  fact  that  the  stiffness  matrix,  [Qij],  for 
each  lamina  is  constant  within  that  layer.  Therefore,  the  stiffness 


matrix  can  be  removed  from  behind  the  integration  but  remains  within 
the  summation  sign.  Also  note  that  mid-plane  strains  and  curvatures 
are  Independent  of  Z  and  can  also  be  removed  from  the  integration  and 
summation.  Thus,  the  force  and  moment  resultants  can  be  written  as: 
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3  3 

®lj  “  y  Z  ^^ijH  ~  ®k-l^ 


The  Aij's  are  called  extensional  stiffnesses,  the  Bij's  are  called 
bending-extensional  coupling  stiffnesses,  and  the  D^j's  are  called 
bending  stiffness.  For  a  linear  analysis  the  mid-surface  strains  and 
curvatures  can  be  expressed  in  terms  of  the  displacements  by 


substituting  Eq.  (8)  into  Eq .  (7)  to  yield 
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(16) 


where  the  non-linear  terms  have  been  left-off. 

STAGS-Cl  Theory 

The  Structural  Analysis  of  General  Shells  (STAGS-Cl)  is  a  computer 
code  developed  by  the  Lockheed  Palo  Alto  Research  Laboratory  to  analyze 
general  shells  under  various  static,  thermal,  and  mechanical  loading 
[25-26].  The  approach  used  by  STAGS-Cl  to  solve  these  problems  is  an 
energy-based  finite  element  analysis.  According  to  the  energy  method,  a 
system’s  total  potential  energy  Is  used  to  derive  Its  equilibrium 
equations  from  which  stability  can  be  determined  by  the  solution  to  an 
eigenvalue  problem. 

A  shell's  total  potential  energy,  V,  Is  equal  to  Its  Internal  strain 
energy  minus  the  product  of  the  external  forces  and  their  respective 
deflections.  This  can  be  written  as  [27); 

^  tKl{<3}  -  d  {ri}  (17) 
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where  Is  the  vector  of  nodal  degrees  of  freedom  of  the  structure, 

[K]  Is  the  structural  stiffness  matrix,  and^Rl^ls  the  structural  applied 
load  vector. 

An  element's  strain  energy,  U,  Is  given  by  [24]: 

1  r  i  (T  [[A]  [B]' 

^  ~7  J  Wo  [B]  [D] 

area  •- 

where  Is  the  mid-plane  strain  and  curvature  vector  given  by  Eq .  (7) 

and  [A],  [B],[D]  are  the  3x3  stiffness  matrices  given  by  Eq's.  (13) 
and  (14). 

In  general,  the  strain  vector  Is  a  function  of  the  mid-surface 

displacements  (u,  v,  and  w) ,  the  first  order  partial  derivatives  of  u, 

V,  and  w  with  respect  to  x  and  y,  and  the  second  order  partial 
derivatives  of  w  with  respect  to  x  and  y.  Bauld  [24]  carried  out  the 
Integration  of  Eq .  (18)  using  the  terms  from  Eq.  (7)  for  found 

that  the  expression  for  strain  energy  Is  comprised  of  three  distinct 
parts.  The  first  part  Is  quadratic,  the  second  part  Is  cubic,  and  the 
third  part  Is  quartlc  In  displacements.  This  result  can  be  written  In 
terms  of  the  appropriate  element  shape  functions  and  nodal  degrees  of 
freedom . 

By  using  a  similar  finite  element  analysis  on  the  element's  external 


{‘l< 


dA 


(18) 


forces  the  element's  external  potential  energy  can  be  obtained.  This 
combined  with  the  strain  energy  gives  the  total  system  potential  energy 
Eq .  (17)  and  can  be  written  In  the  form  given  by  Bauld  [24]  as: 


Where  [I]  and  [J]  represent  nonlinear  stiffness  matrices  in  unknown 
displacements  and  products  of  displacements,  respectively.  For  a  linear 
analysis,  X,  Is  the  proportionality  constant  of  a  convenient  load  level 
used  In  the  equilibrium  Eq's.  (20)  to  solve  for  the  unknown 
displacements.  Also  the  J  matrix,  which  arises  from  the  prebuckling 
rotations,  Is  omitted.  The  quantities  H,  I,  J,  or  equivalently, 

Nlfg,  and  are  calculated  once  based  on  the  equilibrium 

displacements.  Finally,  the  load  proportionality  parameter,  X,  is 
incremented  until  a  sign  change  on  the  left  side  of  Eq.  (21)  occurs, 
signifying  bifurcation. 


(• 
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III.  Finite  Element  Modeling 

This  thesis  will  evaluate  the  stability  of  cylindrical  composite 
panels  that  are  subjected  to  a  simple  shear  load.  Such  a  panel  would  be 
representative  of  an  aircraft  fuselage  skin  panel.  The  STAGS-Cl 
requires  the  user  to  Input  the  type  of  shell  surface  geometry.  To  model 
the  cylindrical  panels  and  to  allow  for  different  radii  of  curvature, 
the  STAGS-Cl  cylindrical  shell  geometry  was  selected.  (See  Figure  3.1). 
In  Figure  3.1  U,  V,  and  W  represent  the  displacements  In  the  X,  Y,  and  Z 
directions,  respectively  and  RU,  RV,  and  RW  represent  the  rotations 
about  the  X,  Y,  and  Z  coordinate  axes,  respectively. 


Figure  3.1  STAGS-Cl  Cylindrical  Shell  Geometry 
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This  study  also  required  the  equivalent  of  a  flat  plate.  To  model 
this,  using  the  cylindrical  shell  geometry  model,  a  study  was  run  to 
determine  What  value  of  R  would  correspond  to  00.  For  a  symmetric 
laminate,  a  radius  of  10,000  In.  gave  values  within  0.02  percent  of 
STAGS-Cl  flat  plate  geometry. 

Boundary  Conditions  and  Loading 

The  boundary  conditions  used  in  the  bifurcation  analysis  were 
selected  to  represent  simple  shear  in  a  flat  plate.  The  boundary 
conditions  (See  Figure  3.1)  assumed  that  the  structure’s  ring  frames  and 
stringers  were  effective  In  restraining  out-of-plane  deflections  (W), 
also  the  top  and  bottom  frames  restrained  the  rotational  movement  about 
the  X-axls  (RU)  and  the  two  sides  stringers  restrained  the  rotational 
movement  about  the  T-axls  (RV).  The  boundary  conditions  are  summarized 
In  Table  1,  where  0  represents  a  fixed  displacement  and  1  represents  a 
free  displacement  along  the  panel's  edges.  These  values  were  used  for 
both  the  pre-  and  post-buckle  boundary  conditions. 

Table  2.1 

Panel  Boundary  Conditions 

U  V  W  RU  RV  RU 

Top  110  011 

Right  Side  110  101 

Bottom  0  0  0  011 


Left  Side 


110 


10  1 


To  Input  a  ainple  shear  load  Into  the  panel,  line  loads  (N^y  and 
Nyx)  vere  applied  to  all  four  edges  (See  Figure  3.2).  They  vere  assumed 
to  act  uniformly  along  all  sides  of  the  panel  which  In  a  flat  plate 
would  represent  a  simple  shear  loading. 


Grid  Size 


The  determination  of  a  finite  element  grid  size  Is  a  trade-off 
between  cost  and  accuracy.  Nelson  [29]  found  that  to  obtain  accurate 
results,  there  should  be  at  least  five  node  points  per  each  half  sine 


wave  of  the  buckled  pattern  In  the  circumferential  direction.  Several 
papera  [2,  3, and  7]  indicate  that  for  the  panel  used  in  this  analysis 
a  grid  of  0.5  z  0.5  in.  elements  represent  a  good  trade-off  between 
accuracy  and  computer  economics.  For  the  12.0  x  12.0  In.  panel,  used  in 
this  analysis,  a  0.5  x  0.5  In.  element  would  correlate  to  a  mesh  of  24  x 
24  elements  or  a  total  of  576  elements  and  625  nodes.  Due  to  the  large 
number  of  computer  runs  needed  to  complete  this  analysis,  it  was  decided 
to  Increase  the  element  size  to  2/3  x  2/3  In.  which  correlates  to  a  mesh 
of  18  X  18  elements  or  a  total  of  324  elements  and  361  nodes.  (See 
Figure  3.2)  Computer  runs  showed  that  this  still  gave  five  nodes  per 
each  half  sine  wave  and  reduced  the  cost  per  run  by  2.3  times.  To 
varlfy  that  this  reduction  would  still  yield  valid  results,  a 
convergence  study  was  run  which  sho%/ed  that  the  difference  between  the 
bifurcation  loads  for  the  324  elements  and  576  elements  is  approximately 
0.24  percent  which  was  well  trlthin  engineering  accuracy. 

Finite  Element  Selection 

The  STAGS-Cl  finite  elements  program,  as  used  in  this  analysis,  uses 
flat  elements  to  model  curved  surfaces.  Figure  3.3  Illustrates  how  a 
cylindrical  surface  can  be  modeled  with  flat  elements.  The  use  of  flat 
elements  presents  conformity  problems  which  are  important  when  dealing 
with  nonlinear  and  stability  problems.  Compatibility  problems  develop 
for  both  the  rotational  and  displacement  degrees  of  freedom  at  the  nodes 

The  rotational  problem  Is  dealt  with  by  assuming  that  the  angle  of 
intersection  a  Is  small,  and,  as  a  consequence,  the  normal  rotation  /32 


Is  Ignored  while  the  conformity  constraint  Is  satisfied  by  letting 

Complete  displacement  compatablllty  along  the  common  boundary 
requires  that: 

cos  (a/2)  -  (w^^^  +  w^^^)  sin  (a/2)  ■  0 

(23) 

(w^^)  -  w^^^)  cos  (  a/2)  +  (v^^^  +  v^^^)  sin  (a/2)  ■  0 

Where  v  and  w  are  the  displacements  In  the  y  and  z  directions, 
respectively.  These  conditions  will  not  be  met  If  v  and  w  are  not 
represented  by  polynomials  of  the  same  order.  To  satisfy  this 
requirement,  additional  degrees  of  freedom  are  added  to  the  element  to 
raise  the  order  of  the  polynomials  representing  the  Inplane  deformation. 
Since  w  Is  represented  by  a  cubic  polynomial,  because  the  strain  energy 
expression  includes  second  order  derivatives  of  the  transverse 
displacement  v.  It  Is  necessary  that  u  and  v  also  be  represented  by 
cubic  polynomials.  This  Is  achieved  by  the  use  of  two  rotations  at  each 
corner  node,  -  v,^  and  u,y,  and  tangential  displacements  at  mid-side 
nodes.  The  difference  between  these  two  rotations  yields  a  shear  strain 
at  each  corner  node,  which  Is  Introduced  as  an  additional  degree  of 
freedom.  Thus,  each  element  has  a  total  of  32  degrees  of  freedom.  In 
STAGS-CI  this  element  Is  referred  to  as  the  5R411.  (See  Figure  3.4) 

A  somewhat  simpler  version  of  the  SH411  element,  the  SH410  element. 
Is  also  Included  for  thin  shell  analysis.  This  element  excludes  the 
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■ld~slde  tangential  dlaplacenents  and  uses  only  an  average  normal  rotation 
at  each  of  the  corner  nodes.  This  restricts  u  to  a  linear  function  In 
the  x-dlrection  and  v  to  a  linear  function  in  the  y-Jirectlon,  Also  the 
shear  strain  is  suppressed  at  the  comer  nodes.  (See  Figure  3.5) 

Because  of  this  the  SH410  element  Is  an  Incompatible  element,  but  it 
still  yields  valid  results  for  some  problems.  For  a  more  thorough 
explanation  of  these  elements  development,  the  reader  should  refer  to 
Reference  {30}  and  [31]. 

To  help  assist  In  the  decision  of  ffhlch  element  to  use,  a  small 
convergence  study  was  conducted  for  both  elements  (See  Figure  3.6  and 
3.7).  The  geometric  radius  tras  selected  to  be  24  In.  for  the  12  In.  x 
12  In.  panel.  Figure  3.6  shoves  the  differences  between  the  eigenvalues 
for  the  SH411  and  SH410  elements  verses  the  number  of  elements  In  the 
grid.  The  SH410  converges  to  a  slightly  higher  value  (about  1.6  percent 
higher  at  324  elements)  than  the  SH411.  The  reason  for  this  Is  that  the 
SH410  produces  a  stlffer  element  than  the  SH411  due  to  It's 
Incompatibility.  Figure  3.7  shows  a  8.2  percent  difference  In  the 
bifurcation  load,  at  324  elements.  The  values  for  N^y  were 

calculated  by  multlpuling  the  eigenvalues  by  the  Internal  resulting 
force  Nxy*  This  force  was  calculated  by  STAGS  based  on  the  applied 
load  and  the  element  selected.  For  this  analysis  a  unit  line  load  N^y 
was  applied  to  the  boundaries.  Therefore,  the  difference  betmen  the 
SH410  and  SH411  values  In  Figure  3.7  was  due  to  the  differences  in  the 
elements  and  how  STAGS  distributed  the  load  through  the  panel.  Since  the 
SH410  suppresses  the  shear  strain  at  the  nodes.  It  requires  a  higher 
shear  load  to  produce  bifurcation. 


R1 furcation  Lodd,  N  (Ib/in) 


R  =  24  in. 
[0/45/-45/90]. 


•  Sti410 
O  SH411 


Even  though  the  SH410  element  does  not  perform  as  mil  as  the  SH411, 
It  was  selected  for  use  In  this  analysis  when  evaluating  the  Influence 


of  hygrothermal  and  curvature  effects.  Part  of  this  reason  Is  cost.  The 
SH411  cost  between  2  and  2.5  times  that  of  SH410  run.  The  other  reason 
Is  that  this  analysis  Is  Interested  In  the  trends,  so  when  the 
bifurcation  load  Is  normalized  with  respect  to  the  room  temperature  and 
zero  moisture  values,  the  differences  In  the  bifurcation  load  reduces  to 
a  ratio  of  the  eigenvalues  and  the  error  divides  out. 


Bifurcation  Analysis  Hethod 


I  (• 


STAGS-Cl  has  two  buckling  analysis  modes.  One  uses  a  linear, 
pre-buckllng  displacements  state  and  Is  referred  to  In  the  STAGS-Cl 
manual  as  the  bifurcation  analysis  with  a  linear  stress  state.  The 
second  method  uses  a  geometric  pre-buckllng,  nonlinear  displacement 


calculation  and  Is  referred  to  as  the  bifurcation  analysis  with  a 
non-linear  displacement  state.  For  this  study  the  bifurcation  analysis 
with  a  linear  stress  state  will  be  used.  This  method  will  calculate 
the  pre-buckllng  displacements  and  rotations,  stress  resultants, 
strains,  and  stresses  as  desired.  It  also  predicts  the  bifurcation 


eigenvalue  and  Che  shape  of  Che  eigenvector.  It  does  not  yield  any 
post-buckling  Information. 


Validation  of  Model 


To  check  the  boimdary  conditions,  loading,  and  how  well  STAGS 


predicts  buckling  with  bending-extenslonal  coupling,  the  finite  element 
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model  ms  compared  to  a  Galerkln  solution  obtained  by  Whitney  [1].  In 
his  paper,  he  presented  an  analytical  solution  for  the  buckling  of 
anisotropic  cylindrical  plates  under  arbitrary  combinations  of  axial 
load.  Internal  pressure,  and  In-plane  shear  load.  He  used  the  Donnell 
equations,  as  they  are  applied  to  laminated  cylindrical  shallow  panels. 
In  conjunction  with  the  Galerkln  method  to  determine  the  critical 
buckling  loads.  Two  of  his  solutions  were  selected  to  check  the  finite 
element  model.  One  related  to  the  effects  of  bending- twisting  coupling 
and  the  second  related  to  the  effects  of  bending-extenslonal  coupling. 

This  comparison  considered  two  12-layer  angle-ply  laminates  with  the 
stacking  geometries  [+453/-453]g  and  [+45^7-45^]  under  pure  in-plane 
shear  loading.  The  panels  were  18  x  18  In.  square  with  a  thickness  of 
h  ■  0.06  In.  In  order  to  keep  the  same  number  of  elements,  the  grid 
spacing  was  increased  to  1  In.  x  1  in.,  which  still  gave  five  nodes  per 
half  sine  wave.  The  buckling  boundary  conditions  also  had  to  be 
changed  In  STAGS  to  the  simple  support  boimdary  conditions  referred  to 
as  BC-1  to  correlate  to  Whitney's  work.  These  boundary  conditions  are: 


At  X  ■  0,  18  in. 

Nn  -  Mn  -  V  -  W  -  0 

and  at  y  ■  0,  18  In. 

Nn-Mn“U-W-0 


(24) 


The  Galerkln  method  satisfies  the  natural  boundary  conditions.  Whitney 
also  used  the  following  ply  properties  for  his  laminates: 


E1/E2  -  lA.O,  Gy2/E2  -  0.5,  1/^2  -  0.25  (25) 

which  are  tjrplcal  properties  of  current  graphite /epoxy  composite 
materials.  These  values  had  to  be  modified  somewhat  for  Input  into  the 
STAGS  program  to  give: 

El  -  14  X  106  psl 
E2  “  i  *  10®  psl 

^  (26) 
Gi2  *  0.5  X  10®  psl 

V^2  ■  0*25 

which  yellds  Whitney's  ratios. 
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Figure  3.8  shows  the  comparison  for  the  model  studying  the  effects 
of  bending- twisting  coupling.  The  two  sets  of  curves  correspond  to  a  + 
or  -  shear  load.  The  symmetric  lay-up  l+453/~453)g  eliminates  all 
from  Eqn’s  (13)  an  (14)  along  with  Aig  and  A2g»  but  leaves  In  the 
bending- twisting  coupling  terms  Djg  and  D2g*  Since  Whitney  used  the 
Donnell  equations  as  shown 
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He  had  to  restrain  the  radius  parameter  In  order  to  keep  the 
assumptions  used  by  the  Donnell  equations  valid.  The  smallest  radius 


Figure  3.8  STAGS-Cl  Finite  Element  vs  Whitney’s  Galerkin,  [453/~A53] 


he  investigated  was  for  an  R  >  60  In.,  which  correlates  to  a  h/R  x  10^ 
of  1.0.  Whitney's  iiork  was  extrapolated  linearly  down  to  a  radius  of 
12  In.  for  comparison  purposes.  This  Is  possible  because  Donnell’s 
equations  are  linear  In  R.  Both  STAGS  elements  SH410  and  SH411 
correlate  to  what  Whitney  found  In  the  range  from  R  ~  60  to  00.  What 
differences  were  found  came  about  because  Whitney  was  able  to  Input  a 
pure  shear  load,  Into  his  equations  and  zero  out  all  remaining 

resultant  forces  and  moments.  With  STAGS  the  Input  of  a  pure  shear 
load  also  generates  some  additional  resultant  forces  and  moments  which 
contribute  to  reduce  the  bifurcation  load.  The  deviation  between  the 
extrapolated  Whitney  and  STAGS  solutions  Is  primarily  due  to  the  use  in 
STAGS  of  the  general  shell  equations: 


The  underlined  terms  are  functions  of  the  radius  which  are  not  found  In 
Donnell's  equation  and  come  into  play  as  the  radius  decreases.  The  use 
of  the  Donnell  equation  for  a  radius  of  12  In.  would  Introduce  an  error 
from  34  -  39  percent  If  compared  to  the  SH410  element  and  from  45  -  54 
percent  If  compared  to  the  SH411  element. 


Figures  3.9,  3.10,  3.11  shov  the  coDparison  for  the  aodel  studying 
the  effects  of  bendlng-extenslonal  coupling  with  Figure  3.11  being  a 
conblnatlon  of  Figures  3.9  and  3.10.  The  unsyametrlc  lay-up 

eliminates  Ajg,  D15.  and  D26  along  with  all  the  B^j's 

except  Bjg  and  B2g,  the  bendlng-extenslonal  coupling  terms,  from  Eqn's 
(13)  and  (14).  One  can  see  that  similar  conclusions  can  be  reached 
from  these  curves.  The  main  difference  Is  that  the  prebuckle 
displacements,  forces,  and  moments  play  a  much  more  Important  role  In 
this  lay-up  no  matter  what  the  curvature  since  they  can  not  be  removed 
from  the  STAGS  runs.  Again,  using  the  Donnell  equation  for  a  radius  of 
12  In.  would  Introduce  an  error  of  around  59  percent  If  compared  to 
the  SH410  element  and  around  79  percent  If  compared  to  the  SH411 
element.  Again,  If  the  correct  value  for  the  bifurcation  load  Is 
needed  the  SH411  should  be  used  as  found  in  the  bendlng-twlsting 
components  and  the  convergence  test. 

The  Influence  of  hygrothermal  effects  can  develop  bending  - 
extenslonal  coupling  In  symmetric  lay-ups.  This  analysis  indicates 
that  the  model  selected  for  this  thesis  should  give  a  good  indication  of 
the  Influence  of  hygrothermal  and  curvature  effects  on  cylindrical 
composite  panels  loaded  In  shear  if  the  curvature  Is  greater  than 
R  -  60  In.  Also  since  Che  bifurcation  load,  N^y,  will  be  normalized  In 
this  thesis  for  analysis  of  the  effects  of  hygrothermal  and  curvature, 
the  use  of  the  SH410  Instead  of  the  SH411  will  not  Influence  the 
outcome  of  the  analysis  when  comparing  a  curvature  down  to  R  >  12  In. 
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Figure  3.9  STAGS-Cl  Finite  Element  (SHAIO)  vs  Whitney's 
Galerkln, 
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Figure  3.10  STAGS-Cl  Finite  Element  (SH411)  vs  Whitney’s 
Galerkin, 


IV.  Evaluation  of  Moisture  and  Temperature  Conditions 


Due  to  their  excellent  performance  characteristics,  advanced 
composite  materials  have  been  gaining  wide  use  in  aerospace  structures. 
However,  It  Is  recognized  that  absorption  of  moisture  and  exposure  to 
thermal  environments  can  have  underslrable  effects  on  the  mechanical 
properties  of  such  materials  [7-21].  Shlrrell,  Halpln,  and  Browning 
[19]  reported  that,  unless  these  effects  are  accounted  for  in  designing 
a  system,  the  service  life  and  reliability  of  a  polymeric  matrix 
composite  may  be  compromised.  Moisture  absorption  affects  the  composite 
In  several  different  ways.  Firstly,  the  resin  swells  causing  a  change 
In  the  residual  stresses  of  the  composite  and  possibly  micro-crack 
formation.  Secondly,  the  resin  may  be  plasticized  thus  causing  an 
increase  in  the  elongation  of  the  resin  near  failure.  This  also  has  an 
effect  on  the  damping  of  the  material.  This  plasticization  is  the 
result  of  the  lowering  of  the  glass  transition  temperature  Tg*  The 
glass  transition  temperature  Is  actually  a  temperature  range  below  which 
the  resin  Is  essentially  brittle  and  above  ifhlch  It  behaves  rubbery. 
Thirdly,  the  Interface  between  fiber  and  resin  may  be  affected  thus 
Influencing  the  composites'  strength  and  toughness.  The  fibers  are  not 
affected  by  either  moisture  absorption  or  moderate  thermal  environments 
which  might  be  encountered  during  a  normal  aircraft  service  life. 

These  changes  In  the  resin  have  been  found  to  result  In  a  decrease  In 
the  tensile  properties  [13]  and  a  reduction  In  the  transverse  and  shear 
moduli  [12  and  15]  of  the  composite  material.  Also  a  slight  Increase  In 
the  longitudinal  elastic  modulus  was  reported  in  Reference  [12]. 


Advanced  resinous  composites  can  absorb  water  through  the  following 


mechanisms: 

1.  the  fiber  -  matrix  interface  (capillary); 

2.  cracks  and  voids  In  the  composite;  and, 

3.  the  resin  (diffusion). 

Of  these  three,  the  primary  mechanism  of  water  penetration  in  large 
well- fabricated  composite  structures  Is  a  rapid  surface  absorption 
followed  by  diffusion  of  the  water  through  the  resin  [19]. 

Prediction  of  Moisture  Absorption 

The  testing  of  composite  materials  to  determine  their  mechanical 
properties  Is  usually  done  at  known  temperatures  and  moisture 
concentrations.  Since  the  primary  mechanism  for  absorption  of  moisture 
Is  through  diffusion.  Pick’s  second  law  of  diffusion  [32]  which,  under 
certain  circumstances  has  been  shown  to  correlate  with  test  data  [12], 
will  be  used  as  the  model  for  diffusion  in  this  analysis.  Pick 
developed  this  equation  In  1855  by  drawing  an  anology  between  heat 
conduction  In  a  solid  and  diffusion  through  a  solid.  Pick's  second  law 
Is: 


dt 


(28) 


where  C  Is  the  concentration  of  the  diffusing  substance  through  the 
thickness  of  the  laminate  as  a  function  of  time  and  distance  through 


the  thickness,  Z  is  the  space  coordinate  measured  normal  to  the 
surface,  K  is  the  diffusion  constant,  and  t  is  time. 


The  solution  of  this  partial  differential  equation,  vith  boundary 
and  Initial  conditions  pertinent  to  the  problem  of  diffusion,  is  shown 
below.  This  is  the  same  solution  used  by  Snead  and  Palazotto  [7]  and  is 
slightly  different  from  that  form  found  in  Section  4.3.3  of  Reference  [32] 
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Where  C  is  as  previously  defined,  Cq  is  the  initial,  uniform  moisture 
concentration  through  the  thickness  of  the  laminate,  Ci  and  C2  are  the 
initial  moisture  concentration  conditions  at  the  inside  (-z)  and  outside 
(+z)  surface  of  the  laminate,  respectively,  and  h  is  the  thickness  of 
the  laminate. 


The  moisture  concentration  distribution  through  the  thickness  can  be 
determined  by  using  this  series  solution  with  a  known  diffusion  constant 
and  prescribed  initial  conditions.  Then  the  reduced  mechanical 
properties  of  each  ply  can  be  obtained  from  appropriate  test  data  by 
assuming  that  the  effective  moisture  concentration  of  each  ply  can  be 


approximated  by  the  calculated  moisture  concentration  at  the  middle  of 
that  ply< 

This  series  solution  Is  a  combination  of  a  steady  state  moisture 
distribution  (the  first  ttro  terms)  and  a  transient  moisture 
distribution  (the  last  two  terms)  which  decreases  with  increasing  time. 
Therefore,  the  accuracy  of  this  series  approximation  Is  dependent  upon 
the  number  of  terms  used  during  the  two  summations.  Snead  [7]  wrote  a 
computer  program  which  was  used  In  this  analysis  to  calculated  the 
solution  to  the  series  approximation.  To  Insure  an  accurate  solution, 
the  program  carries  out  the  summation  until  there  Is  no  change  from  the 
previous  answer.  For  the  CDC  computer,  this  Is  equivalent  to  14 
significant  digits  of  accuracy.  A  version  of  Snead's  computer  program 
Is  shown  in  Appendix  A. 

The  use  of  Pick's  equation  has  certain  limitations  which  must  be 
considered.  The  series  solution  of  Pick's  equation  was  derived  assuming 
a  constant  moisture  diffusion  coefficient  K.  In  reality  the  diffusion 
coefficient  Is  a  function  of  the  laminate's  temperature  and  moisture 
concentration.  The  diffusion  of  moisture  into  a  composite  can  be 
assumed  to  take  place  at  a  constant  temperature  in  simple  cases  since 
the  process  Is  relatively  slow,  requiring  many  months  or  years  before 
the  laminate  achieves  a  state  of  equilibrium  moisture  concentration. 


The  accuracy  of  Pick's  equation  to  model  moisture  diffusion  in 
composites  Is  also  affected  by  rapid  temperature  changes.  The  rapid 
thermal  heating  of  the  laminate  to  temperature  near  the  material's  T 


has  been  found  to  Increase  the  rate  of  oioisture  weight  gain  above  that 
which  is  predicted  by  Pick's  equation  (11,  12,  14,  and  17].  This 
increase  is  believed  to  be  due  to  the  development  of  surface  crazing  and 
cracking  brought  about  by  the  rapid  heating  and  resin  swelling  [12]. 

Pick's  equation  has  been  generally  accepted  as  a  good  initial 
approximation  of  the  moisture  concentration  distribution  through  a 
composite  laminate,  for  simple  cases,  when  the  restrictions  of  no  rapid 
heating,  surface  crazing  or  cracking,  and  assuming  that  a  constant  K  are 
observed  [12,  16,  18,  and  19].  These  are  the  assumptions  being  made  for 
this  analysis. 

AS/3501-5  Mechanical  Properties 

The  material  properties  required  by  STAGS-Cl  as  input  parameters  are 
the  composite's  longitudinal  modulus  Ei>  transverse  modulus  E2,  shear 
modulus  Gi2»  Poisson's  ratio  V2I’  Poisson’s  ratio  l/-^2  relates  the 
strain  in  the  2  direction  to  the  strain  in  the  1  direction  when  stressed 
in  the  1  direction.  Snead  and  Palazotto  [7]  used  the  experimentally  - 
measured  data  for  the  graphlte/epoxy  system,  AS/3501-5,  found  in  Plgure 
8.18  of  Reference  22  to  determine  the  elastic  moduli  as  a  function  of 
temperature  and  moisture  concentration.  The  values  of  E2  and  Gj^2  used 
in  his  work  from  which  intermediate  values  are  linearly  Interpolated, 
are  shown  in  Table  4.1.  Por  a  better  perspective,  Pigure  4.1  sho%fs  this 
data  in  graphical  form.  The  malarial  properties  as  calculated  by  Snead 


and  Palazotto  will  be  used  in  this  work  in  order  to  correlate  the 
find ings. 


TABLE  4.1 


Values  of  Transverse  and  Shear  Moduli 
for  AS/3501-5  [15] 

Transverse  Modulus,  E2 

(psl) 

Moisture 
Concentration 
( percent) 

80  “F 

Temperature 

200 ‘F  250 ‘F 

300  “F 

0.0 

1.41375E06 

1.09475E06 

1.015E06 

1.015E06 

0.050 

1.305E06 

0.9135E06 

0.6235E06 

0.522E06 

1.050 

1.2615E06 

0.841E06 

0.4785E06 

0.290E06 

Shear 

Modulus ,  Gi2  (psl) 

0.0 

0.8555E06 

0.7830E06 

0.6815E06 

0.6525E06 

1.050 

0.8555E06 

0.6597E06 

0.3915E06 

0.1522E06 

The  moisture  and  temperature  Influences  on  the  transverse  and  shear 
moduli  are  clearly  evident  In  the  experimental  data  for  AS/3501-S  shown 
In  Table  4.1  and  Figure  4.1.  The  transverse  modulus  £2  shows  degradation 
both  at  room  temperature  and  elevated  temperatures  while  the  shear 
modulus  Gi2  only  shows  degradation  at  elevated  temperatures.  Also 
vlsable  Is  the  change  In  the  Tg  caused  by  the  absorbed  moisture.  The 
resulting  plasticization  of  the  resin  Is  shown  by  the  Increased 
degradation  In  the  moduli  with  Increasing  moisture  concentration  at  each 
elevated  temperature.  The  longitudinal  modulus  Ej  Is  dominated  by  the 
fiber  stiffnesses  and  hence  is  not  significantly  Influenced  by  changes 
In  moisture  and  temperature  as  are  the  matrix  dominated  £2  and  G]^2 
moduli.  Therefore,  Ei.  Is  assumed  to  be  contant  and  has  a  value  of 
18.85  X  10^  psi  for  AS/3501-5.  Also  for  this  work  the  value  of  V\2 
assumed  to  be  a  constant,  V-^2  “0.3. 

Moisture  Conditions 

The  series  solution,  Eq .  (28),  to  the  Pick’s  equation  requires  three 
moisture  concentration  coefficients  Cq*  Cj^,  and  C2  to  determine  the 
moisture  concentration  through  the  laminate's  thickness  as  a  function  of 
time.  These  moisture  concentrations,  measured  as  a  percentage  of  the 
weight  gained  by  the  composite,  correspond  to  the  initial  moisture 
concentration  In  the  laminate,  the  moisture  concentration  at  the 
Interior  (-Z)  surface,  and  the  moisture  concentration  at  the  exterior 
(+Z)  surface,  respectively,  for  the  cyllndlcal  panel.  Table  4.2  lists 
the  three  moisture  concentration  conditions  considered  by  Snead  and 
Palazotto  [7]. 


Table  4.2 


Snead  and  Palazotto's  Moisture  Conditions 


Cond.  No. 

Co 

Cl 

C2 

1 

0.00 

0.00 

0.0105 

2 

0.00 

0.0105 

0.00 

3 

0.00 

0.0105 

0.0105 

They  used  a  surface  moisture  concentration  of  0.0105  to  correspond  to 
the  material  test  data,  Table  4.1,  that  was  available  for  a  saturation 
moisture  concentration  of  1.05  percent.  (Note:  The  concentration  of 
1.05  percent  relates  to  an  environment  which  has  a  75  percent  relative 
humidity.)  Conditions  1  and  2,  in  Table  4.2,  result  in  an  unsymmetrlc 
degradation  of  the  E2  and  0^2  resulting  In  an  unsymmetrlc 

laminate,  which  %rlll  Introduce  bending-extension  coupling.  Condition  3, 
In  Table  4.2,  is  symmetric  and  will  not  produce  any  bending-extension 
coupling.  Snead  and  Palazotto  found  that  conditions  1  and  2  produced 
similar  results  to  each  other,  therefore,  only  one  un63nBmetrlc  condition 
was  examined  In  this  thesis  along  with  one  symmetric  condition.  Table 
4.3  shows  these  moisture  conditions. 

Table  4.3 

Moisture  Conditions 


Cond .  No . 

CO 

Cl 

C2 

1 

0.00 

0.00 

0.0105 

2 

0.00 

0.0105 

0.0105 
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Along  with  the  coefficients  Co>  and  C2,  In  Eq.  (29),  It  Is 
standard  to  use  a  dlmenslonleas  time  t*,  where  t*  •  K  (In^/sec)  x 
t(sec) /h^(ln)^.  This  eliminates  the  need  to  pick  a  specific  diffusion 
coefficient,  K  and  time.  The  dimensionless  times  used  in  this  analysis 
are  0.0,  0.001,  0.01,  0.1,  and  0.5.  These  values  represent  from  zero 
moisture  absorption  to  a  steady-state  distribution,  respectively.  Tabl 
4.4  shows  how  t*  corresponds  to  real  time  given  a  specific  K  and 
laminate  thickness  h. 

Table  4.4 

Relation  Between  Real  and  Dimensionless  Time 
Real  Time  Real  Time  Dimensionless  Time  t* 


(sec) 

(days) 

0.0 

0.0 

0.0 

3.045E04 

0.35 

0.001 

3.045E05 

3.52 

0.01 

3.045E06 

35.24 

0.1 

1.527E07 

176.24 

0.5 

Note:  These  times  were  calculated  using  K  ■  0.52537E-10  (in^/sec) 
for  an  8-ply,  0.04  thick,  AS/3501-5  laminate. 

The  parametric  equation  used  to  determine  the  value  of  K  ■  0.52537E-10 
2 

(in  /sec)  in  Table  4.4  was 

Kdn^/sec)  -  6.51  exp(-5722/T)  (0.03937)^  (30) 

where  T  Is  the  laminate  temperature  in  degrees  Kelvin.  This  equation 
was  taken  from  Reference  [221.  Figure  4.2  shows  the  moisture 


distribution  through  the  thickness  for  the  five  tlae  values  and 
nolsture  conditions  1  and  2.  This  figure  shows  that  the  panel  will 
absorb  moisture  nonsymmetrlcally,  for  moisture  condition  1.  This  means 
that,  as  time  Increases,  the  panel  will  reach  a  steady  state  condition, 
t*  -  0.5,  where  each  lamina  has  a  different  moisture  content.  The 
symmetric  nature  of  moisture  condition  2  Is  also  shown,  and  as  time 
Increases,  the  panel  will  reach  a  steady  state  condition,  t*  *  0.5, 
where  each  lamina  has  the  same  moisture  content.  The  moisture 
distribution  Is  shown  as  a  continuous  function  In  Figure  4.2,  however, 
as  mentioned  previously,  the  moisture  concentration  Is  calculated  at  the 
center  of  each  ply  and  then  assumed  to  be  constant  through  the  ply 
thickness  In  this  analysis. 


Temperature  Conditions 

For  this  analysis  the  laminate  was  assumed  to  be  at  a  constant 
temperature.  This  Is  possible  since  the  moisture  diffusion  process  Is 
much  slower  than  that  of  thermal  diffusion.  Crank  [32]  Indicated  a 
factor  of  10^  for  the  time  difference.  Four  different  temperatures;  80, 
200,  250  and  300*F;  were  used  In  this  analysis  to  evaluate  the  Influence 
of  a  wide  range  of  temperatures  on  a  composite  panel.  This  range  varies 
from  room  temperature  to  the  operational  survlce  limit.  The  service 
limit  for  graphlte/epoxy  is  approx Imatelly  300"F.  These  temperatures 
were  selected  because  the  material  test  data.  Figure  4.1,  was  taken  at 
these  four  temperatures. 


Laminate  Ply  Orientations 


Two  elght-ply  laminates  will  be  evaluated.  The  two  ply  orientations 
chosen  are  the  [0/45/-45/90]8  and  the  [45/-45]28»  which  correspond  to  two 
of  the  three  orientations  Investigated  by  Snead  and  Palazotto  [7]. 


Re  suits  and  Discussion 


For  the  hygra thermal  Investigation  carried  out  within  this  study,  a 
combination  of  moisture  conditions,  temperatures,  and  times  generates  a 
matrix  of  40  cases  for  each  given  laminate  per  radius.  These  40  cases 
are  broken  Into  two  sets  of  20  cases;  20  cases  for  each  moisture 
condition.  Because  of  the  large  number  of  cases  to  be  evaluated,  all  40 
cases  were  run  for  only  the  radius  of  12  In.  for  both  laminates.  In  the 
analysis  of  curvature,  the  radii  36,  48,  96,  and  10,000  In.  were  run  for 
only  the  [45/-45]2s  laminate,  with  a  limited  number  of  temperature  and 
time  steps.  This  reduced  the  number  of  cases,  but  still  provided 
adequate  Information  to  analyze  the  effects  of  curvature  on  the 
bifurcation  load.  In  this  last  mentioned  analysis,  moisture  condition  1 
was  evaluated  at  300*F  and  time  steps  0.00  and  0.5,  while  moisture 
condition  2  was  evaluated  at  80,  200,  and  300 *F  for  all  five  time  steps. 
A  matrix  of  the  case  numbers  and  corresponding  conditions  Is  shown  In 
Table  5.1  and  a  list  of  the  bifurcation  loads,  Njjy,  for  these  cases  can 
be  found  In  Appendix  B. 

The  computation  of  the  moisture  concentration  distribution  and  the 
reduction  In  the  transverse  and  shear  moduli  iS  done  by  the  computer 
program  listed  In  Appendix  A.  This  computer  program  also  generated  the 
Input  deck  for  the  STAGS-Cl  program,  and  is  a  modification  of  the 
computer  program  written  by  Snead  [7].  It  was  modified  to  allow  for 
shear  loading  and  changes  In  curvature. 


Table  5.1 


Moisture,  Tenperature,  and  Radius  Conditions  Evaluated 


Case 

No. 

Moisture 

Condition 

Temperature 

Radius 

(In.) 

Laminate 

1-20 

1 

a  ,b  ,c  ,d 

12 

[0/45/-45/901 

21-40 

1 

a,b,c,d 

12 

[45/-45128 

41-60 

2 

a.b,c,d 

12 

t0/45/-45/90] 

61-80 

2 

a  ,b,c,d 

12 

I45/-45]28 

81-91 

2 

a  ,b  ,c  ,d 

24 

[45/-45l2s 

92-93 

1 

d 

24 

[45/-45128 

300-305 

2 

a  ,d 

36 

f45/-45J28 

100-110 

2 

a,b,d 

48 

t45/-45]2s 

111-112 

1 

d 

48 

t45/-45]2s 

121-128 

2 

a,b,d 

96 

[45/-45]2s 

131-132 

1 

d 

96 

I45/-45]2s 

141-147 

2 

a,b,d 

10,000 

t45/-45]28 

Note:  The  cases  with  a  R  ■  12  In.  were  run  for  all  five  times  and  four 
temperatures,  however,  the  remaining  cases  were  run  at  a  mixture 
of  time  steps  and  temperatures  depending  on  the  moisture 
condition  and  ra-llus. 

Notation: 


Moisture 

Condition 

1  —  Co  “  Cj  -  0.00 

C2  “  0.0105 

2  —  Co  -  0.0 

Ti  “Co  *  0.0105 


Temperature 

(•f) 

a  —  80.0 

b  --  200.0 
c  --  250.0 


Reduction  In  Bifurcation  Load 


The  evaluation  of  the  data  for  a  panel  loaded  In  ahear  Indicates 
that  the  degradation  of  the  E2  ^12  moisture  and 

temperature  effects  results  In  a  reduction  of  the  panel's  bifurcation 
load,  Nxy  Figures  5.1  and  5.3  show  the  results  of  the  STAGS-Cl  runs 
for  the  two  12  In.  by  12  In.  panels,  with  a  radius  of  12  In.,  while 
Figures  5.2  and  5.4  show  the  results  found  by  Snead  and  Palazotto  [7] 
for  a  similar  panel  acting  under  compressive  loads.  In  these  plots 
Nxyorjg  *^orlg  bifurcation  loads  for  a  panel  at  80*F 

and  a  dimensionless  time  of  t*  ■  0.00  when  loaded  In  shear  and  axial 
compression,  respectively.  These  values  are  unaffected  by  either 
temperature  or  moisture  degradations.  The  actual  shear  buckling  values 
are  shown  In  tabular  form  along  irtth  Individual  plots  for  each  time  and 
temperature  series  In  Appendix  B. 

As  was  expected,  the  panel's  bifurcation  load  decreased  with 
Increasing  temperature  and  absorbed  moisture.  Comparing  Figure  5.1  to 
Figure  5.2  shows  a  similarity  in  the  trends  for  the  reduced  bifurcation 
loads  for  the  [0/45/-45/90]s  lanilnate  considering  shear  and 
compression.  A  point  should  be  made  that  even  though  the  percent 
reductions  are  similar  the  actual  values  for  the  bifurcation  loads. 
Table  5.2,  differ  considerably,  l.e.,  for  room  temperature,  zero 
moisture  Nx  *  514.8  lb/ In  compared  to  ••  123.4  lb/ In  for  the 
[0/45/-45/90]8  laminate’  difference  In  Figures  5.1  and  5.2  Is  that 

the  curves  for  the  shear  loading  are  shifted  down  below  those  for  axial 


Reduced  Bifurcation  Load 


•  sH 


2 

Dimensionless  Time  (Kt/h  ) 


Figure  5.1  Degradation  In  Nyy  For  the  l0/45/-45/90]g  Laminate 


Reduced  Bifurcation 


Reduced  Bifurcation  Load 


Figure  5.3  Degradation  in  K^y  for  the  [45/-45]2s  Laninate 


Reduced  Bifurcation  Load,  N  /N 


Table  5.2 


Comparison  of  Bifurcation  Loads  (lb/ in) 
at  80“F,  t*  -  0.00,  and  R  -  12  in. 

Axial 

Laminate  Compression  [7]  Shear 

l0/45/-45/90)s  514.8  123.4 

[45/-45]2g  ^28.9  160.9 

loading.  For  example  at  a  temperature  of  SOO^F  and  moisture  coniitlon 
2,  the  reduction  in  shear  was  26  percent  while  for  axial  compression 
it  was  21  percent,  a  difference  of  5  percent.  A  similar  analogue  can  be 
drawn  from  Figures  5.3  and  5.4  for  the  [45/-45]s  laminate,  except  that 
the  curves  for  the  shear  loadings  are  shifted  up  above  those  for  an 
axial  load.  For  example  at  a  T  •  300 *F  and  moisture  condition  2,  the 
reduction  for  shear  was  35  percent  while  for  axial  compression  it  was  43 
percent  a  difference  of  8  percent.  To  illustrate  this  better.  Figures 
5.5  thru  5.8  show  these  comparisons  for  each  moisture  condition  and 
laminate.  Figures  5.5  thru  5.8  also  show  that  the  laminates  perform  as 
designed.  That  is,  the  rate  of  reduction  for  the  ■I45/-45]2s  laminate 
loaded  in  shear  was  less  than  that  for  the  [0/45/-45/90 ]g  laminate 
loaded  in  shear  %»hlle  the  rate  of  reduction  for  the  l0/45/-45/90]g 
laminate  loaded  in  compression  was  less  than  that  for  the  (45/-45]2s 
laminate  loaded  in  compression  for  a  given  temperature  and  moisture 
conditions . 

As  with  Snead  and  Palazotto's  [7]  work,  the  reductions  in  the 
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Figure  5.8  Axial  vs  Shear  Loading  for  Moisture  Condition  2, 

(45/-45I2S 

58 


bifurcation  loads  are  significant  considering  that  diaenslonal  changes 
due  to  resin  swelling  were  not  Included.  The  maximum  reduction  In 
for  each  laminate  and  moisture  condition  Is  summarized  in  Table  5.3  for 
the  radius  of  12  in. 

Table  5.3 

Percent  Reduction  in  Bifurcation  Load, 
at  300*F,  t*  -  0.5,  and  R  -  12  In. 

Laminate  Moisture  Condition 

1  2 

[0/45/-45/90]g  16.7  25.6 

[+45/-45]2s  17.8  34.5 

Another  similarity  to  Snead  and  Palazotto*s  [7]  work  Is  that  for 
moslture  condition  1  the  reduction  in  the  bifurcation  load,  %y>  was  not 
as  great  as  It  was  for  moisture  condition  2  even  though  moisture 
condition  1  causes  the  initially  symmetric  laminate  to  become 
unsymmetrlc,  which  introduces  bending  -  extension  coupling.  The 
symmetric  moisture  condition  2  has  a  much  greater  Influence  on  both 
laminates.  This  is  due  to  an  over  all  general  reduction  In  the  material 
properties  as  moisture  is  absorbed  into  the  panel  symmetrically. 

Effects  of  Curvature 

The  I45/-45]2s  laminate  was  examined  at  six  different  radii  (12, 

24,  36,  48,  96,  and  10,000  in.),  considering  moisture  condition  2,  to 
evaluate  the  effects  of  curvature  on  the  bifurcation  load.  This 
laminate  was  selected  because.  In  comparison  to  the  [0/45/-45/90]g 


laminate.  It  exhlblta  a  greater  reduction  In  the  reduced  bifurcation 
load  for  a  given  moisture  content  and  temperature.  Figure  5.9  shows 
the  general  overall  effects  of  curvature  on  the  bifurcation  load.  That 
Is,  as  the  panel's  radius  Increases  the  bifurcation  load,  N^y 
decreases.  Also  shown  In  this  figure  are  the  effects  of  temperature  and 
moisture  absorption  on  the  bifurcation  load.  Though  the  only 
temperature  shown  Is  300*F  with  times  absorption  set  at  t*  *  0.0  and 
0.5,  it  can  be  stated  that  for  any  given  radius  the  Increase  in 
temperature  and  moisture  content  reduces  the  bifurcation  load.  A  better 
Illustration  of  the  effects  of  curvature  Is  shorn  in  Figure  5.10.  This 
figure  shows  that  curvature  does  influence  the  percent  reduction  In  the 
bifurcation  load.  This  is  different  from  Snead  and  Palazotto's  [7] 
findings.  They  Investigated  three  radii,  12,  24,  and  48  In.  and  came  to 
the  conclusion  that  for  a  panel  under  axial  loading  the  reduction  In  the 
panels'  bifurcation  load,  did  not  significantly  vary  from  those  for 
the  12  In.  radius  panel.  They  therefore  concluded  that  the  results 
obtained  for  the  12  In.  radius  panel  should  be  valid  for  any  radius. 
Figure  5.11  shows  how  this  conclusion  can  be  reached  by  only  examining 
panels  of  radii  12,  24,  and  48  in.  If  the  radii  96  and  10,000  in.  are 
added  to  their  work,  the  effects  of  curvature  become  apparent.  The 
overall  effect  of  curvature  Is  to  decrease  the  influence  of  temperature 
and  moisture  cn  the  panels  bifurcation  load.  This  means  that  using  the 
reduced  bifurcation  load  results  for  the  12  In.  radius,  or  any  other 
radius,  as  a  valid  approximation  for  another  radius  will  yield  Invalid 
results.  Thus  In  general,  the  effects  of  curvature  must  be  accounted 
for.  One  needs  to  compare  results  for  the  reduced  bifurcation  load  for 
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Reduced  Bifurcation 


Reduced  Bifurcation  Load 


Che  same  radius  of  curvature.  Note  chat  although  only  a  few  points  are 
plotted  In  Figure  3.11  the  trend  Is  still  quite  visible. 

Looking  at  Figures  5.10  and  5.12,  one  can  see  that  the  Influence  of 
curvature  Is  also  a  function  of  Che  laminate's  temperature.  These 
figures  are  plots  of  the  reduced  bifurcation  load  for  the  [45/-45]2g 
laminate  for  moisture  condition  2.  Figure  5.10  has  the  plots  for 
temperatures  80  and  300®F,  while  Figure  5.12  Is  a  plot  for  the 
temperature  of  200®F.  It  is  apperent  from  these  plots  that  for  a 
temperature  of  80®F  the  effects  of  curvature  can  be  neglected,  and  Chat 
as  the  temperature  Increases  to  300 ®F  the  effects  of  curvature  becomes 
more  Important  and  cannot  be  neglected.  As  the  temperature  Increases, 
the  effects  of  curvature  tends  to  decrease  the  Influence  of  moisture 
and  temperature  on  the  reduced  bifurcation  load.  This  difference  can 
be  seen  best  by  comparing  the  maximum  reduction  In  the  bifurcation 
loads  for  the  panels  of  radius  12  In.  and  10,000  In.  The  12  In.  panel 
has  a  35  percent  reduction  In  Che  bifurcation  load  while  the  10,000  In. 
panel  has  only  a  12  percent  reduction.  This  Is  a  23  percent  decrease  In 
the  effects  of  moisture  and  temperature  on  the  reduced  bifurcation 
load.  This  reduction  Is  also  apparent  at  Intermediate  temperatures, 
(Figure  5.12)  but  to  a  lesser  extent.  Examination  of  the  same  panel 
but  at  moisture  condition  1  also  shows  the  effect  of  curvature  on  the 
reduced  bifurcation  load  (See  Figure  5.13).  However,  the  overall 
effect  Is  not  as  drastic  as  It  Is  for  moisture  condition  2.  Note  that 
this  plot  Is  for  the  worst  case  only,  which  Is  at  a  temperature  of 


Reduced  Bifurcation  Load 


An  interesting  anomaly  In  Figure  5.10  Is  that  the  decrease  In  the 
reduced  bifurcation  load  Is  not  sequential  In  Increasing  radii.  This 
can  be  seen  by  comparing  the  curves  for  24,  36,  and  48  In.  at  300*P. 

The  panel  with  a  24  in.  radius  falls  above  the  panels  with  a  36  In.  and 
48  In.  radius.  This  Is  contrary  to  what  the  theory  indicates  In  Eq. 
(16),  the  strain-displacement  relationships,  where  the  curvature  shows 
up  as  a  linear  dependent  variable,  which  Indicates  that  the  curves 
should  be  sequential.  Since  Figure  5.10  is  a  plot  for  a  normalized 
bifurcation  load,  ^y^^xyorig*  question  arose,  "Is  this  anomaly 
actually  in  the  bifurcation  load  or  is  It  due  to  the  normalizing 
factor?”  Looking  at  the  bifurcation  load  first.  Figure  5.14  shows  that 
as  the  radius  increases  the  bifurcation  load  decreases  for  a  given 
temperature  and  moisture  condition.  This  was  also  shown  In  Figure  5.9. 
Both  of  these  figures  show  that  the  bifurcation  load  Is  following  the 
trend  Indicated  by  the  theory.  However,  Figure  5.14  also  shows  that 
the  rate  of  decrease  Is  dependent  on  the  moisture  content  for  the  given 
temperature.  To  correlate  Figure  5.14  to  Figure  5.10,  the  bifurcation 
loads  were  translated  to  a  value  of  150  lb/ In  so  that  each  curve  would 
have  the  same  maximum  (See  Figure  5.15).  This  figure  clearly  shows 
that  the  plot  of  the  bifurcation  loads  does  not  show  the  anomaly  shown 
In  Figure  5.10.  Therefore,  this  anomaly  must  be  due  to  the  normalizing 
factor. 

All  the  plots  for  the  reduced  bifurcation  loads  were  normalized  to 
the  room  temperature,  80*F  and  zero  moisture  concentration  bifurcation 
load,  **xyQj.j^g*  bifurcation  load  Is  dependent.  In  general,  on  the 

prebuckled  bending  moments  and  membrane  forces,  with  the  curvature  of 


the  panel  governing  which  load  la  more  dominant.  For  a  flat  panel  the 
dominant  loads,  which  catise  buckling,  are  the  membrane  forces  and  as  the 
panels  curvature  Increases  the  bending  moments  become  Important.  What 
Is  being  seen  in  Figure  5.10  Is  a  transition  In  the  dominant  cause  for 
buckling  from  membrane  to  bending.  This  Is  vlslable  since  the 
normalizing  factor  Is  the  relative  shear  bifurcation  load. 

Prebuckled  Displacement ,  w  and  Eigenvector  Charac terlstlcs 

Figures  5.16,  5.17,  and  5.18  show  that  there  Is  a  similarity  In  the 
appearance  of  the  prebuckled  displacements,  w,  when  comparing  the 
effects  of  temperature  and  moisture  for  a  panel  with  a  12  In.  radius. 

The  primary  difference  being  one  of  deflection  magnitude.  The  terms 
’Max.'  and  'Min.'  In  these  figures  correspond  to  the  maximum  and  minimum 
out-of-plane  displacement  caused  by  the  unit  line  load  applied  on 
the  panel's  boundaries  while  the  term  'Contour  Step  Size'  Indicates  the 
Increment  between  contour  lines.  As  the  radius  Increases  there  is  a 
change  In  the  panel's  prebuckled  displacement's,  w,  pattern.  The  area 
of  negative  displacement  becomes  greater,  for  Increasing  radii,  due  to 
decreasing  effects  of  the  prebuckle  bending  moments. 

The  eigenvectors.  Figures  5.19  thru  5.21,  are  also  not  dependent  on 
the  moisture  conditions  and  temperature,  but  they  are  dependent  on  the 
curvature,  similar  to  the  prebuckle  displacement.  The  eigenvector  for 
the  panel  %rlth  a  12  In.  radius  exhibits  five  half  sine  waves  while  a 
10,000  In.  radius  exhibits  two  half  sine  waves.  This  also  Indicates 
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Figure  5.16  Contour  Plots  for  the  Prebuckled  Displacement,  w,  [45/-45J2 
at  80*F,  t*  -  0.0 


that  the  energy  needed  to  buckle  the  panel,  with  a  12  in.  radius,  is 
greater  than  that  for  a  10,000  in.  radius.  Additional  displacement, 
and  eigenvector  contour  plots  for  representative  cases  of  moisture, 
temperature,  and  curvature  are  included  in  Appendix  C. 


VI  Conclusions 


The  following  conclusions  can  be  made  for  cylindrical, 
graphite/ epoxy  composite  panels  with  various  radii  of  curvature  subject 
to  moisture  exposure  and  elevated  temperatures,  when  loaded  In  simple 
shear . 

1.  The  results  for  the  STAGS-Cl  finite  element  analysis,  using  flat 
plate  elements,  compare  well  to  solutions  obtained  by  Whltney[l],  using 
the  Galerkln  method  with  the  Donnell  strain  displacement  relations,  for 
relatively  flat  shells*  (Radius/thlcknest  ^  1000) 

2.  The  bifurcation  load  of  a  composite  panel,  with  a  resin  material 
whose  elastic  moduli  are  reduced  by  absorbed  moisture  and  elevated 
temperature,  will  degrade  %ri.th  Increasing  moisture  concentrations  and 
temperatures. 

3.  The  trend  for  the  reduction  In  bifurcation  load  of  a  composite 
panel  subjected  to  a  simple  shear  load  is  comparable  to  that  found  for 
an  axial  compression  load  at  a  given  radius. 

4.  The  extent  of  the  degradation  In  the  bifurcation  load  Is 
Influenced  by  the  moisture  concentration,  the  temperature,  the  panel's 
ply  orientation,  and  panel's  curvature.  At  300*F,  a  radius  of  12  in., 
and  a  symmetric  moisture  %relght  gain  of  1.05  percent,  the  I0/45/-45/90]s 
panel  experienced  a  25.6  percent  degradation  and  the  [45/-45]2s  panel 
experienced  a  34.5  percent  degradation. 


5.  The  cylindrical  panel's  bifurcation  load  is  influenced  by  the 
ply  orientation,  and  curvature. 

6.  The  bending  -  extension  coupling  induced  by  the  unsymmetric 
initial  moisture  condition  did  not  significantly  influence  the 
bifurcation  load. 

7.  Increasing  the  cylindrical  panel's  radius  decreased  the  panel's 
bifurcation  load  and  can  significantly  change  the  moisture-  and 
temperature-induced  degradation  characteristics  for  a  given  temperature 

8.  If  curvature  is  an  Important  variable  in  a  given  problem,  the 
use  of  the  normalized  bifurcation  load,  with  respect  to  the  membrane 
forces,  can  lead  to  conclusions  that  are  invalid  if  one  is  trying  to 
model  the  rate  of  change  in  the  bifurcation  load. 


The  main  program  calculates  the  moisture  concentration  at  the 
mld-polnt  of  each  lamina  using  Pick's  Second  Law  of  Diffusion.  The  main 
program  calls  the  following  subroutines. 

Subroutine  Header  —  Creates  tape6  which  contains  the  matrices  A,  B 

and  D  and  the  Qij’s. 

Subroutine  Calcl  —  Calculates  the  reduced  transverse  and  shear 

moduli  and  the  value  for  V21  for  AS/3501-5, 
given  the  moisture  concentration  in  the  ply. 

Subroutine  Calc2  —  Calculates  the  values  for  Qjj  ,  Ajj,  ®od 

Dij. 

Subroutine  Stagsl 

and  —  Creates  tape7  the  Input  for  STAGS-Cl. 

Subroutine  Stags2 


Appendix  A 


Computer  Program 


PACCItAP  MClC  (INPUT  tOUTPUl*T«F(Ct1«PE1) 


C  AkTHONT  0.  STHANi  CIE-S^O 

C  TPIS  PRCGRAN  IS  A  NOOIFICATICA  OF  CAC  OCVEIOPEO  BT  JAPES  N. 

C  SAEAO. 
t 

r 

C  THIS  PRC6RAH  CALCULATES  THE  POISTURE  CONCEATRAT ION  THROUGH  THE 
C  THICKNESS  CF  A  COFPOSITE  LAHINATE.  THIS  PRECRAH  ALSO  CALCALATES  THE 
C  PLT  STIFFNESSES  AND  UVARIAM  PRCPEFTIESt  THE  FLT'S  REDUCEA 
C  STIFFNESSES,  AND  THE  LAMINATE'S  CA',  CB!*  ANC  ED!. 

C  USING  THE  CALCULATED  REDUCED  HOOULI*  INPUT  FILES  FCR  THE  SJAGS 
C  FINITE  ELEMENT  DUCKLING  LQAC  PRC6RAM  ARE  ARITTEN  FCR  A  PAAEL 
C  LOADED  IN  SIHFLE  SHEAR. 

C 

C  REF.  'THE  HATHENATICS  OF  DIFFUSION*  BT  JOHN  CRANK*  SECiNO  EDITICN* 
C  CLAREOON  PRESS*  OFCRD*  ISTE. 

C  REF.  ■INTRODUCTION  TO  CCFPCSITE  PATERIALS*  BY  STEPHEN  A.  TSAI  AND 

C  H.  THCHAS  HAHN*  TECHNCFIC  PUBLISHING  CO.*  1980. 

C  REF.  "HCCHONICS  OP  COMPOSITE  MATERIALS*  BT  ROBERT  M.  JiNES* 

C  MCORAN-HILL  BCON  CtMPANT*  I97E. 

C 

C  BASIC  MOISTURE  DIFFLSICN  COUATICN  -  PICK  EOUATION  <J.  iRANK) 

C 

C  K«S.O.  OF  C  URT  2  :  F.C*  OF  E  URT  T 

C  nhere: 

C  K  s  MOISTURE  DIFFUSION  COEFFICIENT 

C  S.O.  =  SECOND  PARITAL  CERIVATITE 

C  F.C.  =  FIRST  PARTIAL  CERIVATINE 

C  NRT  s  UlTH  RESPECT  TO 

C  C  =  SPECIFIC  POISTURE  CONCENTRATION  IN  LANINAt 

C  7  =  SPACIAl  CCORCINATE  THROUGH  LAMINATE  THICMESS 

C  T  =  TIME 

C 
C 

I  notes: 

C  1.  THIS  PROGRAM  IS  SET  UP  FOR  ASA3S01  GRAPHITE  EPOXT  lUING 
C  MATERIAL  PROPERTIES  FCRM  TSAI'S  TERT 

C  2.  the  temperatures  are  in  UNITS  CF  DEGREES  KELVIN. 

C  3.  the  temperature  cistributicn  through  the  thickness  js  assumec 

C  TO  BE  CONSTANT  AND  EQUAL  TO  ONE  OF  THE  FOUR  SPECIFItO 

C  TEMPERATURES  DEFINED  AS  ElTEMPfAM. 

C  A.  TAPEG  CONTAINS  THE  OUTPUT  FILE  ANC  TAPET  CONTAINS  T*E  INPUT 
C  FILE  FOR  THE  STACS  PROGRAM. 

REAL  K 

COMRON/MAT/2(20I*T(  I  0  I *TNCNCIM( ]0I *ELTEMP ( A ) *L *NO*H*K *iL*LLL *TE MP * 
•CO*C1*C2*RAO*PI*TCMP2(2  0I*C(A,S*8 )*E2(A  *5*6>*U2I(  A*S*84* 
•G1?(A*!*RI«THCTA(8) *011 (A*5«ei*612(A,S*8)* 

‘  «02?(A  tS*8  l•066(A*S*FI*l'l<4*S*(l* 

•U2CA*5*RI*U3fA*S*F» ,UA(A«S*PI*uS«A*«*e)* 

•0R11(A***8I*0R12( A*S*8l*0RlG(A*T*8>*aR22(A*B*B}* 
*0R2€<A,S*8I*QRGG(A,S*8  )*C1*ICASE*HLAM, 

•All(A*SI,A12(A*SI*AlG|A*;i*A2?(A*S) *A !( ( A * E ) * AGG( A  * SI  * 
•Rii(A*si*Bi2(A*;i*Ric(A*!i*e:2<A*:i*6;((A*si*B6e(A*s>  * 
■01I(A*SI*D12(A*S)*DlG(A**l*C:2(A*SI*0:t(A*SI*DGG(A*S> 

C  THESE  ARE  THE  FOUR  TEMPERATURES  AT  NHICH  TEST  DATA  FOR 
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I*. 


i 


I 


'V 

i' 


C  »S/^5P1  IS  «VI1L*BL(  IN  IS*!*!  TtHl 
ELIE"PU»  =  3C0. 

ELTCNP«2li;€6. 

ELIENPCJIsJS*. 

ELTtHP(^»=*22. 

C  *S/3501  Cl  POCUII 
ElslB.e'EOt 
C  7CPO  CUT  THE  »BR*TS 
DO  1  lll=lt« 

DC  1  1=1. E 

DC  1  tlsl.e 

2<Ll>=THET*fLLt=0.0 

01  lILLLtL  .111=012 (ILL «L«tt>  =  G22<LLL*L*ll>  =  066fLUiL«LLJ=0.0 
UMILL.L.LL)  =  U2(LLI  .L.LII  =  U3(ILL  .L. LL  >=U4( LLL.L.LL  >=U5^LLL «L  tLl  > 
•=0.C 

0R11CLLL«L.LLI=0B12<LLL*L .LLI  =  CP1C.(LLI  .ltLL»=0.0 
aB22ILll*L«Lll=eR  26  ILIL.L.LL  >=0R<6(I.LL*l*LL)=(.t 

1  CGNTINUE 

00  4  1=1.10 
«  TCtl  =  TI«CNCIP<ll=0.0 

PBIBT  «.  IC«SC.Cft.C1.C2.H.HLAP.BC.K 
C  INPUT  eCClBlNC  CASE  NC.tlEt 

PRINT  •«  •  INPUT  CASE  AC.  C1S»  ■ 

REAC  •.ICASE 

C  INPUT  CC.C1.C2  IN  UMTS  OF  PCRCEM/100  fEX.-l.CS  =  0.011 
C  Cl  IS  THE  CONCENTRATION  ON  THE  INSICE  OF  THE  SMEtl 
C  C2  IS  THE  CONCENTRATICN  CN  THE  CUTSICE  OF  THE  SHELL 
C  CP  IS  THE  ORIGINAL  CONCENTRATICN  IN  THE  LAPINATE 
PRINT  INPUT  C0.C1.C2  • 

REAC  •.C0.C1»C2 

C  INPUT  LANINAE  THICKNESS  AND  NLPBER  CF  PLIES 

PRINT  «.  •  INPUT  LANINAE  THICAAESS  ANC  NO.  CF  PLIES  • 

REAC  •.  MLAP.NO 
C  CALCULATE  lANlNATE  THICKNESS 
H=Nt»HLAK 

C  CALCULATE  THE  CENTER  OF  EACH  lAAINAE 
7(1  )  =  HLAN/2.0 
00  2  L=2.NC 

2  2(L)=2<L-1I«HLAP 

C  INPUT  THE  PLT  OPIENTATICNS  FRCP  INSICE  TO  tUTSlOE 

PRINT  •,  •  INPUT  PLT  CRIENTATICNS  FRCP  INSIDE  TO  OUTSIIE  • 

OC  :  1=1. AC 
READ  •.THETACLI 
S  CONTINUE 

C  INPUT  THE  SHELL  RADIUS  IN  INCHES 

PRINT  •.•INPUT  THE  SHELL  RADIUS  lA;  INCHES* 

REAC  •.  RAD 

C  DIFFUSION  lENPERATURE  =  300  K  (80  F» 

TtPP=300. 

C  REF.  TSAI  FOR  EQUATION  8. AO  TC  CALCLLATE  K  FCR  GRAPHITE/EPAKT 
C  N=R0*E«P(-EC/RT I  (8. AO) 

C  K=6.51*EKP(-5T22/TEAP(DEG  KELVIN))  (FIG  B.T) 

K=6.51«EKP<-5T22FTE AP)«C.03S3T««2 

C  INPUT  NC.  OF  tine  INTENVALS  AT  MHICF  MOISTLRE  DISTNIBUTION  IN  THE 
C  LAMINATE  HILL  BE  CALCULATED 

PRINT  «•  •  INPUT  NO.  CF  TINE  INTERVALS  ■ 

REAC  *.  RTIPE 

C  TRPUT  NCNOIPENSIONNL  TIPE  INTERVALS  (K*T/H**2I 
PRINT*.  •  INPUT  TINE  INTERVALS  (P»T/H**2>  * 

DC  3  Ltl.NTlHE 

3  REAC  •.INONDIP(L) 

P1=3.1A15B2€53«BBTB; 

ICASE=ICASE-1 


' 


f' 

4. 


L 

t- 

C'-; 


V'.- 

u- 


I' 
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flO-RlO  878  THE  BUCKLINO  OF  CONPOSITE  CVLINORICHL  PMELS 

CONSIDERING  ENVIRONNENTflL  EFFECTS(U>  AIR  FORCE  INST  OF 
TECH  URIQHT-PflTTERSON  AFB  OH  SCHOOL  OF  ENGI. . 
UNCLASSIFIED  A  D  STRAH  DEC  8S  AFIT/GRE/'RR/'85D-15  F/G  28/11 


f ' 

r'. 

END 

MICROCOPY  RESOLUTION  TEST  CHART 

N*’iONA|  RnRf  au  Of  t^lANfaRPS  A 


00  lOlO  LLI<1«4 
00  1000  LsltNTIOC 
T(ll=TkONOIf>(LI*H**£/R 
ic«sc=ic/isr«i 
coil  moto 

C«Ll  ST«6SI 
00  3000  Ll^lfRO 
TC)*P2(LL>  =  CLTtMP(LLl) 

C(lLL«l«Ll)=0. 

IFIKLI.CQ.Q.I  EO  TC  12 
CfllLtltLLI=Cl«tC2-Cl)*2f ll»/H 
CTFPPlrO. 

CT€PP2=0. 

M  =  0 

11  0=0«1 

IFIinS(-K*N**2*Pl**2<T<L>/H**2l.fi£.*C0.l  EC  TC  12 

CTCPPl=(2/Pn*(C2*CCS(h*Pl>>Cl»/P*SIkCR*Pl*ZILL>/HI 

CTCPP1=CTEPP1*CPPI-P*R**2*PI**2«1{0/P**2> 

C(LlL«LiLLf=C<LlL«L  .LLMCTEPPl 

P=<P-n/2 

CTt'*P2rl««CC/Pl  >/(2*P«l  t*<!R((?«P«l  )*P  1*2(11 1 /PI 
CTrPP2  =  CTEPP2*l.lP(-p*(2*p«li**:*Pl**2*1(LI/H**2) 
C(llLtl*LL)=C(LLLtL  tLll«C1fPP2 
N:k«l 

CTCPP3r(2/pn*(C2*CCS(P«PI)-Cll/P*Stll(«*Pt*2<LLI/H) 
CTEPP3sCTEPP3*EpP<-«*H**2*Pl**S«T(LI/P**2> 
C(LLL«L«LL)=C(LLLtLtLll«CTEHP3 
GO  TO  11 
12  CONTINIE 

IF(C(LlL«LtLL).LT«0.0>  C( ILL *t *Ll 1=0. 0 
Z2i2(Lll/H 

CALL  CALCl(C(LLL*L*LLI«TePPZ(LL l«C2(LLl «LtLLI »U21 (ILL »i»LLI » 

•  G12(LLL«L*LLI I 

102  F0PPAT(TlC«ll*T30f*tI>l  OEFIAES  PATERIAL  NC«*I 
VAITC(7«102>LL 

100  F0RPAT(4K,tlt2PtFE.4«2l*F6*«»l«*F7.5«lkiC11.5illtF8.5*JS*Cll.S*lP« 

•  I3ilX«I3i 

101  FOP*AT(T3t*lA*8SC06t*tF10.7»***tF10.1***.0*6*l.«««F10.4«*tl* 

•  «l-2«l 

URITE (7«101IU21(LLL  «L«Lll iC12(LLl tL tLL I tC2 (LLLtL «LL I 
WRITC(EflOClLLtZ2tZ(LL)tC(LLLfLtLLI tC 3 (LLLtLtLl I t U21 ( LlLt L tLL I t 
•G13(LLLtLfLLItNtP 
30CO  CONTIAVE 

CALL  CALC2 

CALL  stags: 

1000  CONTINUE 
STCP 
ENC 


SUPNCUTINE  KAOEN 
REAL  R 

COPPON/PAT/?(?C)«T(  101  (TNCRDIP(10l«ELTEPP(AliL«NO«HtR«i.LtLLL«TEPP« 
•C0tClfC2«RAD(PI«TCPP2(:ritC(A«*  tPI*C2(4«StPliL':i(4*5*BJt 
•G12(«»St8 ItTMETAIPI t011(AfSt8l*C12(«i*«8l* 
•022(««EtS)t0E£(4tStP)tLl(4fS»eit 
•U2((«*tRltU3(4«!fPtfU4(4*S*PI«tl(A*lt8)t 
•0Rll(4tS(8ltaR12(4tEf8l«QRl€(4t;*8lt0R:2(4«5«Blt 
•0R2((4<3*ei«eR6E(Af 1 18 l«EltICASttHLAR< 

•All(A,SltA12( AtSI tAlE(*tSI*A:2(A«!l«A:E(A|*>«AE(( At5l t 
•'Bll(4»«l»P12(Atll*816(A(SI»e22(A«:itB:E(A«:i«8EE(A«S)  t 
•011(4,S>,012(4.;i«01E(4t«l«C:2(4««lfD:((4tt)iOA((4t5) 
lor  F0RPAT(1N1*//I 
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10 1 

•T53t*G12»tT€»»«l»«»TI8»»P»«/»lll»*«lfc.  l•#132^•»COUCtO•  t 

•  T»«t»«£CUCE0"«T;3.»SCDl;CEC*./l 
10?  rORP«T<//l 

10«  ro"P*TIT3»*C»SE  NC.  •|IE»*  SHE«»  l»lST»eilITT  IN  COHPOSilE  P»NELJ"/ 
•T;,«C«ICUL«TI0N  of  FOISIUFE  E1S1PIBUT1CN  IN  "tiai 

••PIT  L*PIN4TE*/TSt«FLT  C«tCNT*TICN  < • ,f 5 .1 »•» • t F5 . 1 » • #2 tF 5. 1 »• • •• 
•F?.lt",»iF5  .l,«t»«F? .1, •,•,£;.! »•»•.£ i.l #•)••/ 

*T5*’N0ISTUPE  CIFFtSlON  CCEFFICIFFT  s  •# 

•Ell.at/iT;.«OIFFUSI<N  TENPERFTOEs  •.r7.2**  CNtL*IN»»» 

•/tTE»»PlT  THICKFESJ  =  C.005  1 N . • » /»T5 *»C0  =  •tFO.Si*  Ci  =  *• 
•f9.5#»  C2  =  ••Fa.5t/««  ORICINAL  NOCULI-« F»T£ ••£!  =  18.«5E0«»» 

•  •  C2=1.A1?T5C06  021  =  0.02250  €  12=0 .851*E06  <PSI»*»/»T54 

••HONOIPENSICNAL  TIFE  =  •.£ 11 .5 »5» ••! 1 FE  ISEO  =  ■tEll.lt/l 

10  HRITCfCtlOO) 

HR  ME (Elio A  IICASEfNCiTHETA(lli1PEtA(2>iTHETA(31 iTHETAC  4lt 
•THETA(«)|TPETA(E) iTFETA(7liTHEtA(8l|K|lEFPtC0iCliC2i 

•  TNCNOIFfLI  iKLI 
?0  CONTINOE 

URITECEilOl  I 

RETIRA 

ENC 

SUBRCUTINE  CALCK ( 1 7EPP ^£2  iU21 iE 12 > 

DIPENSICN  S(Ai3)iSS(Ai2lf«(2> 

NsO 

IF(7EPP.EO.?OO.I  A:l 
IFdEPP. £0.366.)  k=2 
IF<1EPP.E0.30A.)  F=? 

1F(TEPP.E0.A22.)  6=4 

1F(N.EO.OI  PRINT  *i*  ERROR  IN  CALC* 

C  REF.  TSAI  FOR  AS/3S01  C2  AKO  612  PCCULl 
S(lill:1.41?75C06 
S(ll2):l.?05E06 
S<lf3)=1.2615E06 
S(2il)=l.e54T*C06 
S(2i2)=.01?5E06 
S(2i3)=.841E0£ 

S(3il)=1.015E06 

S(3i2):.6:35EC6 

S(3f3l=.47e5E06 

S(4ill  =  1.01?C(l6 

S(4f2)=. 522(06 

S(4  f3)=.250E0E 

SS(lfl)=.055SCO6 

SSI1.2I=.8555E0E 

SS(2il)=.7e3CE0E 

SS(2.2)=.E507«Ce6 

SS<?.l)  =  .6f  15E0( 

SS(3f2l=.351SE06 

SS<4il l=.ES25E0E 

SS(4i2)=.15225C06 

«(1)=0.00 

)((2)=0.005 

R(?l=0.0105 

012=0.300 

IF(C.6E.C.0C5I  €0  TC  17 
SLCPC=(S(Ri2)-S(Nil  I)/<»<2I-»( I )) 

B=S(N.l)-SLCPr*»<l) 

E2=SL0PE«C*P 
021=  U12*E2/18.eSE0E 
60  TO  20 
10  -CONTINIE 


SLOPC  :  (S(ll«3l-S(Ht2n/(ll(!»«)it3l> 

B=S(llt2)-SLCPC*SI3l 

t2=StOFt*C«e 

U31:Ui;*C2/ie.BStF( 

20  COkTINCC 

SL0PC:<SS(M«2t-SSIK«llW(ll<2>-litl>> 

B:SSIII«l)-SLOPt*l(l  ) 

C13:SL0PE*C«e 

PCTURN 

CNO 

SUEPOUTINC  ST«SS1 
PC«L  R 

COPPON/H«T/2(20>«1<  lot  «1hCROIPOO>*eL1EPPC4)*LiROiH«RtiLtLLL  tIEPP 
•CO«CltC2«P*0«PI*TePF2f20)*CC4*;«ei*E2l4»5tet*L21f4tSiSJt 
•612 (««■«« It THCT*< at  toil l4t0taitC12(4(>«8)t 
•022l4t!ieit066MtStPltllf4t;t8tt 
•U2(4t3t8  ttU3(4 t^tPI tU4(4t:t(ltVE<4t:tCI t 
•oaii(«fSiaitQai2Mt;tei  tOPiEiat^tSi  toa;2C4  tStOit 

•  0R2C(4i*te)t0R6E(4t‘ta  ttEltlCASEtPL *81 

•All  («tElt«12(4t:i**16(4t!lt«22t4»!) »t2f (4»:)t«€(l«tS> t 
•Bll(«t:itai2(4t3l tai6<«tSttP22(4t!l«B2{(4t;ite(C(*tS) t 
•Oll(4tEltO12(«t0l tC16(4t!ttO:2l4t;i »D:{(«t!>tO(fC«tSI 
10c  FOaa«T<*/JCP*t/*4CSTa«btCP24000C.  C«S121S*STa«ktEMtS47wl*t 
•/•/USER*! /•CMRPEEt* ."t 

•  /■RfL  =  240000.-t/*SE1TLf  700  l.■t/*«TTPCPtST•ES^/UI>:D8200iO.■t 
•/•«TT«CHtST«GS2/UR=r820CB0.*t/*STA6Sl**t/*RETl1(RtST«€Si**t 
•/•REL=240000.»t/*SlACS2,**/ 

•*RE1UR»tST«GS2.*t/ 

••SURITEtT«PE21tC«SE*tl4,-P.«Cl:Cri0472.*t/ 

•*SURITCtT«Pt22tCASC*tl4.*St«Cl:CC104T:.*t/ 

••/ECR«I 

BRl1C<1tlOO)lC*SE«ltOOf 1C«SE«1C10 

101  F0RR«T(T3**C4SE  NO.  •«IEt*  SHCRR  INSTtCILITT  IN  CONPOS.jTe  PANEIE* 

•  /•C  PIT  ORIENTRTIOR  < ■ tPS.l »*«*«ES.l t 't ‘t FS.l t*t * tFS. Jt* t” t 
•P5.1#«#»tr5.1t“t" tFE.lt*t*tF;.l t"»*t/ 

••c  BCINDAPT  CCROITICR  -  TCP  -  ItVtPV.RW  FREE  •»/ 

••C  5ICES  -  UtV.PltRH  FREE  •?/ 

••C  ECTTOP  .  RV.Pk  FREE  •%/ 

•■C  PRREt  RADIUS  =  •.C4.3.»  IRCFES"t/ 

••C  ELEPENT  TTPE  =  41C"./ 

•  ■C  TRANSVERSE  AND  SHEAR  PCOILI  PEOt'CEC  DUE  TC  POISTURl  ARC  TEPF* 
•F"C  ORIGINAL  MODULI  FOR  AS/3S01  6R APHITE/EPOIT  IS*tF 

•"C  Cl=18.PSE0b  E2:1.4137SCtG  U21=.022S  612:C.8SS5ElC  fPSlI*./ 

••C  MOISTURE  OIFFUSIOR  COEFFICIENT  =  •.Ell. 5."  AT  ■.Fb.l**  0E6 

•  RCIVIR«./»C  CO  =  ■tEll.Et"  Cl  =  •tCll.St*  C2  =  ■tCll.*./ 

•■C  Cl  IS  INSIDE  SHELL!  C2  IS  CLTSICE  SHELL"./ 

••C  NONCIPENSIONAL  UPt  =  •.Ell.*.’  lIPE  ESECt  =  ’.Eli.S./ 

••C  LAMINATE  TEMPERATURE  =  ’.FE.!.’  IRElVINt’t 
HRITEET.lOl IICASE.TFCTA(l>tTHETAC2I.THETAI3>*THETAI4l. IHETACSt. 
•THETAIEl.THETAETl.TFETACa ltRAO.M.TEMP«CC.Cl*C2.TNONOIMAl) .TCLI. 
•ELTEMPfLLlI 

102  FORMAT(T10."l.l.ltO.O*.T30t’tB>l  LINEAR  EIFURCATION  AN^L"./ 
•TlO.’l’.T  J0.’$n-2  1  SHELL  UMT  ISEO"./ 

•TIO.’O.O.l.C’.TJO.’lfl-’  8  MATERIALS!  !  1  SHELL  4ALI  PRtPERTT’./ 
•T10.’l.’.T3C."tC-l  STARTING  lOAC  FACTOR"./ 

•TlO.’l tO.CSO- .T30 t"SO>2  NO.  CF  E IGENVAIUES I IMAI  HO.  CPi  SECS"./ 

•  Tn."l".T;C."*0-3  EICENVAIUE  PER  CLISTER  TC  BE  COMPUTE.,"./ 
•T10."14.19".T30."SF-1  19  MCWS.  15  COLUPNS") 


uRfTC(7(ie;i 


ERC 

SUeRCUTINf  ST«6S2 
REOI.  *t 

COPI'OI«/R«7/2(?OI|1t  10)tTNCNCIPIie>»CLlErPt<)tLtltOtNtRtlL*Lll«TCrPt 
•C0*Cl*C2«P«DtPl*TCPF2(2C)iC<4«‘*8l«C2(«i9t8>itZll«tS*ej« 

•612(4  ei  «011(4  t5te>tC12(4»:»B  li 

•022(4  «:*B  )ta6C(4,5,*)«lll(4tEtPI« 
•U2(4*Stfl)«ll3(4t;*PltU4(4,«,(l*l*(4«:*(l* 

•aRii(4(:tP»«0Pi2(4t!«a  >tOPif(4«!«a>«eP22(4»:«8>* 

•  eR2((4<S«8>«0R6((4t!«8  >tC]«IC«!f tHL*P» 

•411(4|«)t412(4.*l •( U(4*;>t*:2(4««|,«2C(4*E)t*66(4*S>t 
•Rll(4i;)«R12(4**l«P16(4«*t*e:2(4««t*BS((4«;) tB6((4t9>* 
•Dll(4t^)*C12(4«*l«01£(4tS>«C22(4«;t«C2((4*E»tC66(4«S> 

105  F0RP*T(T10#"ltl»*iIJ«T25t»»R-l  1  0*11  COEF.  NO.*  SEN  L*»ER  0*tli* 
•12#"  PLIES«»FT10t*l  •  .0C5»  •»F(.2,»  ll(-2"» 

MRITE(Ttl05»N0t*0.TFET*(l> 

N0TEPP=N0-1 

•N6n=12*lF0/(R*0^PI ) 

DO  1  J:2tN0iePP 
1  URlTEdtlOE  IJ«THET*(JI 
lot  F0RP*T((9,12«"»  .005*  •»F<.2I 
107  F0RN«T(T9*I2«*«  .OO!*  •«F6.2«*  *R-2*> 

HRlTE(7*107INr*TH£Tf(h(l 

10*  FCRP»T(T10.«5-»T30.»*P-1  CTLIRC*1C*L  SPELL  CEORETRT'*/ 

•  T5.-0.*  12. «  0.«  •♦H0.7.".*»ES.;»«  »R-2»"*/ 

•T10»"l»»7I0,"*P-5  SPELL  0»LL  CCEf.  *0.  1*  0*7*  1»BLE»*4 
•710»"410«»T30*»*N-1  ELCPER1  17FE**/ 

•710*«0*0»0«n"«T!0*«IP-l  e0LNC«F7  CONS7RAIN7S  OCFINCO  0*  NE»7  C*»0» 
•»/T10*«ll0»eil»«T3e,«»P-2  ICF  -  l,O.RO*RWsFREE»  0 *ROsF.i»ED" »/ 
•710««J J0*lCl»»T3C»«7P-2  4ICP1  SICE  -  L*V»Rll*ROsFREE*»/ 
•T10.»000*Cll"*730*"IP-2  BCIICP  •  RV*RO  =  FREE**/ 
•71flt"110»101«*T30*«IP-2  LEFT  SICE  -  U»»tRU»RO  s  FREE**# 
•T10**1«*T30***6-1  RC.  CF  LC*C  STSTEPS**/ 

•T10*«l*4*0«*T30*«*0-2"*/ 

•T9*«-l.*2*2*l*0«*T3C*«*0-3  TCF  LINE  l(»C"*2 
•T10*-1.*3*1*0.1S"*T;0.-*0-3  RItPT  LIRE  LO*D**/ 
•TlO*«l.*2*S*19*(l»*TIO*-10-3  PC110N  LIRE  LC»D»*/ 
•T9*»-1.*I*1*0*1"*T3C««»6-;  LEFT  LINE  ICFD"*/ 

•T10*»l*l*l*C  •*7I0*«*R-1  CllPtT  CORTRDl*./ 

••/ECR"» 

MRITCdtlOei  •RCB*PFO 

RETIRR 

ENC 

SUBROUTINE  C*LC2 

C  TPIS  SUBROUTINE  C*LCUL*TES  THE  PLT  **0  L*FIR*TE  STlFFNESSCi 
RE*1  R 

OIPENSICN  TTHET*(2a I *22  (201 

COPPON/H«T22(20I«T( 1 0  I  * TRCNCIP( 10  I *EL TEPPI 4 ) *L*RO*H*K*i.L*LLL* TE PP 
•C0.Cl*C2*R*D*PI*TEPF2(2OI*C(4*5,PI*t2(4*5.P».U2I(4*5*8J* 

•612(4  *5 *n  I  .THE T*( 01 *011 (4  *5 *0  I •( 12( 4 *5  *0  I* 

•022(4*5*01*066(4*5*01 *11(4*5*01* 

•U2( 4 *5*8 )*U3( 4*5 *P)*U4( 4 *5*F  l*L5(4*5*0)* 

•ORl 1(4 *5* 8 )*0R12( 4,5*8  I *eRlt(4*E*ei *0822(4*5*81* 
•0R2t(4,5*RI*0R6((4*5*e  I  (E 1  * IC*5 E *HL*H * 

•«11(4*5I**17(4*5I** It(4*5l,*22(4*51*«2((4*5t**6E(4*5)* 
•Bll(4*5>*B12(4*5l*B16(4*51*e:2(4*51«B:((4*5l*B6E(4*SI  * 

•Oil (4*5 1  *C1 2(4*51*0  I6(4*5l*t;2(4*5l*0It(4*51*Ctt(4*51 


C  U12  FOR  AS/liri  (RRPHI1E/EFCRT 
Ut2:|.3O0 
00  1 

C  C«LCUL«TIRG  FLT  STIFFRESSES 

OllELLl  fL*JI=El/ll.-U12*U21(LLl  «L«J>» 
ei2«LLl*ttJI=U12*E2(LLL«l tJ)F( 1 .*01 2*U21 IlLl *1* J>  I 
a22<LLl«L«JI=C2fLLL  il • J 1/ (l.«U12*U21 III t *1 (J >> 

1  Q6E(LLltl.tJ)sGi;(LLliL«<lt 

C  CRLCUERTIRC  TFC  NEOUCCC  FL>  STJFFFESSCS  ISINC  THE  INFRRIAFT 
C  FROPERTIEE  APFROtCH  OF  TS*I  RRO  f«C«RO  (PCF.  JONES) 

00  2  JsltNC 

UlClLLiL«J)=(3.«all(LLEtL*J>«3*«e22CLLl«LtJT*2.*012(LLl«L<J)* 

•4.*0(EfLLl«L«Un/a. 

U2(lLL«L««l)=(011(LlltL«JI-022(LLl«LtJ>>/2. 
U3(LLlfLiJ>=(ail<LLl«lt<fl«022(tll  tL t J l•2. «012(LLL tLt  J >*4«* 
•e6EILLltLfJ)>/P. 

U4CLLLtL*J)=(all(UltLtJ)4a22<LLl  tL*J)«£«*012ILLLtLtJ)*4.* 
•a6EiLii*L(jn/e. 

2  US(LLLiLtJ>  =  t011(LLl tit J ><022  lilt  tL *J  )>2.*B12 IL IL *L * J ) *«.• 
•aSECELltLf JI)/0. 

C  CHANGING  THE  SIGN  OF  THETA  TO  BE  CCASISTART  WITH  THE  TSAI/FACANO 
C  FCRHULATICk  AS  SHOWN  IF  JCFCS 
DO  3  J:l«NC 

3  TrHETA(J»=>THCTA(J>*3.1AlS926S441/180. 

C  CALCULATING  THE  REOUCEC  STIFFNESSES 
C 

00  4  j:I*NC 

0RlI(lLl«L*JI:Ul<LLl«L4j)«U2(LLltLtJ>4CCS(2.*TTHCTA(J)J«U3tUL« 
•Lf  J|4CCS<4«*TTHETA(wn 

0R12(lll«LtJ)3U4fLLl»L*J)-U3fLll*L»J>4CCS(4.*TTHETAtJ)J 

0R22tllL«LtJT3Ul<lLl«ltJT-U2<Llt«l«J)4CCSf2.*TTHETAIJ)J«U3(LLL* 

•L«J)*C0S(4.«TTHETA<J»| 

ORlGEllLtLf J>3-«5*U2fLLl»l«J»*SIN<2.*TTHETACJ>T<U3ILLL»l*J>* 
•SINI4.*TTHE1A(JI) 

«R2E(Lll«L*Jt3-«5*u!(LLL*l«J)4SIN(2**TTHETA(J)>*U3(LLL«lt J>* 
•SIN(4.*TTHETA<JM 

4  OPGSILLlfL  tJUL'SILLl  «L*J)-U3fLLl«Lt  J)*CCSC4.*TTMET  AEJIJ 
C 

C  CALCULATING  THE  Z  COORCINATE  PER  FIG.  4-1  IN  JCNES 
C 

Z7<lT=-NO*HLAF/2. 

H0TENP=H0*1 
00  «  J=2tN0TEPP 
;  zz(j)=zz(j-i)«HiAp 

c 

C  CALCULATING  THE  (Alt  Cf!f  ANC  tC! 

C 

All <LLL  fL)  =  A12<LLLtL)  =  A16<LLL«L  T:A22<LLL«L  UA 2E CLLL«L i : A66CLLL « L I 

•=0.0 

DO  E  J=ltNC 
ZTEPP=22«J*1»-ZZ«J) 

AliaLLfL)  =  All(LLLfl)«ORlllLlLtt.J>*ZTEPP 
A12(LLL«LI=A12(LLL«l)*gP12<LLLtL«J»*ZTCPP 
AlGCLLl*LI=A16fLLltLI*QPlE(Lll«l*JI*ZTtHP 
A22ILLL«L»3A22(LLLtL  »*0R22(lLl «L t J»*ZTEPP 
A26ILLL  iL  1= A2&ILLLiL)*GR2Cf LLLtl«Jt«Z1LPF  ' 
A66aLL«Ll3A66(LLLtL)*0RG6(LlL  »L  tJ>*ZTENP 
E  CONTINUE 

Bll(LLLtL»3ei2(LLl«ll=BlGILlL*L»3B22ILLLtL>3B26ILLLtL):e6GCLLL.LI 


DO  1  J:ltkC 

flll(LLLtL>  =  f>ll(LlL«LI«QRtl<LLL«lt>l>*71CPP*.: 
ni  ?ILlltL  >  =  ILL  tl  >*0)11 2<LLL«l*«)l«71tPP*.‘ 
qt€(LLL«LI=ete<LLL«l >«0RtC<LLL«l*<l)*2TCPP*.: 
P22(Lll«LI=P22(LLLfLI«aP22(LtL«l*d>*21EPP*.: 
B.26(LLltLt=e26(Ul«ll«aP26<Lll«l*J)*21tPP*.5 
B66(LLL  «L>  =  e6efLLLiL l«e»6£iLLL il td) •21CPF* .£ 

I  CONTINIC 

011(LLL«LI-C12(LLL*l >:C16(lLL»l  1=02211 ll »L > =026CL LLtL I -06&( LLL « I > 
•=0  .0 

DO  e  d=ltNC 

2TCPP=22(d*l>*«3-22CdM«3 

Dll<LLl«LI  =  Cll(LLL*ll«0i)ll(LLL*l*d)*2TtBP/3> 
D12CLLL«LI=O12aLLtlt*0<>12(LLltL*d)*2TEPP/3> 
D16(LLLtLI=Cie(LLL«ll*ORl(<LLL*ltd)*23EnP/3> 
D22»tLltL»=C22ULt«ll*0«»22«Llt»l.dl«2TtPP/3. 

02eiLLLfL>  =  C2£(LLL*l>*0P2«(LlL  tl *d>*2TCPP/3. 

D6e(LLLiL l=06t(LLl«ll*«R66(Lll«l*dl*23EPP/3* 
e  CONTINUE 

IOC  FORP»T<//.T6»«PLT«»T10»«Oll«.T2;t"Ol2".T3R.«O22*tT«bt«Af6"» 
•T5Rt»TNET*"tT70.«TE*F“/» 

WRITCtE  «ior  I 

101  FORP/IT(T6«I2«T10tCll.Etl22«Ell.Etl3RiEll.StTREtEll.StT£e* 
•ril.S«TT0«F7.2l 

DO  0  d=l«NO 

0  URITCie  tlOl  KltOll  aiL«l*«ll«012(lLl*l«..  >tC2riLLL  tLfdltOtllLLL  tL  tdlt 
•THET*(d»t7EPr2<*» 

102  FORP»T(//.T«.,»PlT«»T10t*Ul"»722.»O2».T3R.*U3«.TR6t«UR»#73e.»UE»/) 
URITC(til02t 

DO  10  d=ltNC 

1C  URITEieuOl  idfUl(LLl*Lfdl«U2(Lll«L*d»«ES(LlLtL«d)*U*(Lil»L*d)»UEtL 
•Lt  »l *dl 

(0  10!  f0RP«T<lBl//tU««Pl7»,T10.»«Bll*»T22»*CP12«»T3*t*«*l«%T06»»«R22»» 

•TS8.«O*2«"»T70««OBtF»/» 

WRnC(E«10!> 

104  FORP«TIT6tT2«T10«ei l.C*T22tCll.*«T34«Cll.!*T4(>Ell»S»Tte*Ell.£« 
•T73,C11,5I 

DO  11  d=ltNC 

II  WRITC((flO«IJiORll(llLfltdltOR12(ll.LtLtdl*ORiefLLL«L*dJtBR22ILLl*L 

•  t>ll«aR2EfLLltL*d)tOR6€<LLltl«dl 

105  FORP*T<//,Tfc»»*ll=  •tE11.5tT30t"Bll=  • t E 1 1.5tT54»"DH =  •tEll.St/ 

•  T6t»»12=  •»E11.5»  13tlt"B12=  "tE  11.5tT54  i«012=  •»Ell.5»/ 

•T6t"*l6s  ■«E11.5»730««E16=  • fF 1 1 .5» 154 t"Clt=  •»E11.5t/ 

•T6f»*J2=  •fE11.5,T23.»e22=  • .E 1 1 .5, 754 .•C22=  •.Ell. 5*/ 

•T6f«»2f=  •tE11.5»T3C.»e26=  • »E 1 1 .5t 15 4  ,■026=  ■,E11.5,/ 

•76, •466=  •,F11 .5, T30,*fi66=  •,E 1 1 . 5, 754 , *066=  •,E11.5» 

URI1C<6,105U11(LLI,LI  ,P11<LLL,1),D11  (lll,L>, 

•  ai2ail,LI,B12(LLL,ll,C12CLll,l  >, 

•  416fLLL,LI,B16(lLL,l),C16(llL,l  I, 

•  a2  2{LLl,LI,e22<Ll.(.,ll,C22(lll,l  >, 

•«26(LLL,LI,P26<lLl,ll,D26fLLl,L  I, 

•466 CLLl,L I *B66(LLL,l  l,066(LLL,l  ) 

RETIRN 

ENC 
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Summary  of  STAGS-Cl  Runs 


80.0 

123.3802 

123.15182 

122.51071 

121.93600 

121.58738 

1. 00000 

.9982 

.99295 

.98829 

.9864 

200.0 

116.44961 

116.05624 

114.87964 

113.69288 

113.10469 

. 94386 

.94064 

.93109 

.92148 

.91672 

250.0 

113.20656 

112.57167 

109.80117 

107.42864 

105.83039 

.91933 

.91240 

.89167 

.87071 

.86026 

300.0 

112.97977 

111.88204 

108.44757 

104.72933 

102.39487 

.91571 

.90681 

.87897 

.84883 

.83315 

Case  No. 

21-40  (45/-45]2s  ^  * 

12.0  in. 

Laminate 

Nondlmenslonal  Time  (Kt/h^) 

Temperature 

(Deg.  F) 

0.00 

0.001 

0.01 

0.1 

0.5 

160.85196 

1.00000 

160.68808 

.99898 

160.23854 

.99619 

159.84078 

.99379 

159.68968 

.99277 

154.14556 

.95831 

153.828690 

.95634 

152.88125 

.95045 

151.38482 

.94114 

150.40302 

.93511 

148.69028 

.92439 

147.99595 

.92008 

145.89194 

.90699 

142.20926 

.88410 

139.64275 

.86815 

147.30035 

146  .  36451 

143.33539 

137.08092 

132.76267 

Bifurcation  Loads 


Case  No.  41-60  I0/45/-45/90)b  R  -  12.0  in. 


I  (a 


Laminate 

Temperature 


Nondimensional  Time  (Kt/h^) 


(Deg.  F) 

0.00 

0.001 

0.01 

0.1 

0.5 

Nxy/*^xyori8 

80.0 

123.38024 

122.97095 

121.79026 

120.87806 

120.61371 

1.00000 

. 99668 

.98711 

.97972 

.97758 

200.0 

116.44962 

115.62477 

113.09808 

110.77044 

109.90825 

.94383 

.93714 

.91666 

.89780 

. 89081 

250.0 

113.20656 

111.61128 

106.05130 

100.93193 

99.04338 

.91933 

.90461 

.85955 

.81806 

.80275 

300.0 

112.97977 

110.54358 

102.73265 

94.94916 

91.77837 

.91570 

.89596 

.83265 

.76957 

.74387 

Case  No. 

61-80  (45/-45]2s  R  “ 

12.0  in. 

Laminate 

Nondimensional  Time  (Kt/h^) 

Temperature 

(Deg.  F) 

0.00 

0.001 

0.01 

0.1 

0.5 

*^xy^^*yorig 

80.0 

160.85196 

160.42377 

159.69079 

159.10577 

158.92862 

1.00000 

.99828 

.99382 

.99008 

.98898 

200.0 

154.14556 

153.27435 

150  .  94  923 

147.48502 

145.60686 

.95831 

.95379 

.93933 

.91777 

.90608 

250.0 

148.69028 

146.89746 

141.47279 

132.37836 

126.76072 

.92439 

.91411 

.88036 

.82376 

.78881 

300.0 

147.30035 

144.73534 

136.19195 

118.99368 

105.48780 

.91575 

.90066 

.84749 

.73977 

.65519 

90 


Bifurcation  Load 


Case  No.  81-91  I45/-45J28  ^  “  24.0  In. 

Laminate  Nondlmenslonal  Time  (Kt/h^) 

Temperature 

(Deg.  F)  0.00  0.001  0.01  0.1  0.5 

_ 

80.0  103.33830  103.12701  102.57902  -  102.07056 


1. 00000  .99796  .99265  -  .98773 

200.0  99.36737  -  -  -  .96157 

.96157  -  -  -  .91910 

300.0  96.12327  94.59905  89.63421  81.26078  73.17607 

.93018  .91543  .86739  .78636  .70812 


Case  No.  92-93  145/-45128  *  “  24.0  In. 

Laminate  Nondlmenslonal  Time  (Kt/h^) 

Temperature 

(Deg.  P)  0.00  0.001  0.01  0.1  0.5 

”f/^*yorlg _ 

80.0  103.33830  (from  case  81) 

1.00000 


300.0  96.12327  -  -  -  88.27289 

.93018  -  -  -  .85421 


Case  No.  300-305  I45/-45J28  *  “  ^6.0  In. 

Laminate  Nondlmenslonal  Time  (Kt/h^) 

Temperature 

(Deg.  F)  0.00  0.001  0.01  0.1  0.5 

_ 

80.0  82.45344  -  -  -  - 

1.00000  -  -  -  - 


300.0 


76.47433  75.32610  71.41175  63.21319  55.46706 


Bifurcation  Load 


Laminate 
Temperature 
(Deg.  F)  0.00 


Case  No.  100-110  (45/-45128  ^  ”  ^8.0  In. 

Nondlmenslonal  Time  (Kt/h^) 


0.001 


”^y(”*yorig_ 


67.68884  67.55834  67.21995 

1.00000  .99807  .99307 


66.89306 

.98824 


I  L* 


200.0 


300.0 


65.12717 

.96216 


62.88856  61.94464  58.77294  52.30344 


.92908 


.91514 


.86828 


.77270 


Laminate 
Temperature 
(Deg.  F)  0.00 


Case  No.  111-112  I45/-4532s  ^  *  ^8*0 

Nondlmenslonal  Time  (Kt/h^) 


0.001 


Nxy/^3 


*yorlg_ 


62.13212 

.91791 

46.22196 

.68286 


80.0  67.68884  (from  case  100) 

1.00000 


300.0  62.88856 

.92908 


Case  No.  121-128  t45/-45J2s  ^  “  ^6.0  In. 


Laminate 
Temperature 
(Deg.  F)  0.00 


Nondlmenslonal  Time  (Kt/h^) 


0.001  O.Ol 


”*y^^fyorig_ 


56.88047 

.84032 


200.0 


300.0 


45.87038 

1.00000 

44.11947 

.%183 


42.71049  42.04119 

.93111  .91652 


-  42.14792 

-  .91885 

39.88750  36.29267  32.82597 

.86957  .79120  .71562 
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Bifurcation  Load 


Case  No.  131-132  I45/-A5]28  *  *  9^*0 

Laainate  Nondisensional  Tine  (Kt/h^) 

Temperature 

(Deg.  F)  0.00  0.001  0.01  0.1  0.5 

_ 

80.0  45.87038  (from  case  121) 

1 . 00000 

300.0  42.71049  -  -  -  38.84975 

.93111  -  -  -  .84695 

Case  No.  141-147  I45/-45]28  ^  “  10,000.0  in. 

Laminate  Nondinensional  Time  (Kt/h^) 

Temperature 

(Deg.  F)  0.00  0.001  0.01  0.1  0.5 

5[y^5jy2Eii _ 

80.0  26.89203  -  -  -  - 

1.00000  -  -  -  - 

200.0  26.04990  -  -  -  25.34861 

.96869  -  -  -  .94261 


300.0 


25.78407  25.43476  24.50785 


23.68764 


Reduced  Bifurcation  Load 


Reduced  Bifurcation  Load 


APPENDIX  C 


Prebuckled  Displacement,  w  and 
Eigenvector  Contour  Plots 


The  terms  'Max.'  and  'Min.'  in  this  appendix  correspond  to  the 
maxlnvm  and  minimum  out-of-plane  displacement.  For  the  prebuckled 
displacement,  w,  this  out-of-plane  displacement  was  caused  by  a  unit 
line  load,  N](y>  applied  at  the  panel's  boundaries.  The  term  'Step' 
indicates  the  Increments  between  the  contour  lines. 


Prebuckle  Displacements ,w  Eigenvectors,  w 

Figure  C-1  Case  1,  16,  and  20  Contour  Plots  {0/45/-45/*)0 ) 


Prebuckle  Displacements ,w  Eigenvectors,  w 

Figure  C-2  Case  41,  56,  and  60  Contour  Plots  {0/45/-45/90]g 


Prebuckle  Displacements,w 


Eigenvectors,  w 


Figure  C-6  Case  300,  301,  and  305  Contour  Plots  [45/-45)2 


Prebuckle  Displacements ,w 


Eigenvectors,  w 


Figure  C-9  Case  151,  144,  and  147  Contour  Plots  (45/-45]2 
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Abstract 

' ' 

An  analytical  study  was  conducted  to  determine  the  influence  of 
moisture,  temperature,  and  curvature  on  the  bifurcation  load  of 
cylindrical,  composite  panels  subjected  to  a  simple  shear  loading.  Two 
laminate  ply  orientations,  ( 0/45/-45/90]g  and  145/-45)^,  were  analyzed 
for  six  radii,  four  temperatures,  and  two  initial  moisture  conditions. 
■The  eight-ply  composite  panels  were  assumed  to  be  manufactured  from  a 
.graphite/epoxy,  AS/3501-5.  To  evaluate  the  influence  of  moisture  and 
temperature,  the  transverse  modulus,  E^,  and  shear  modulus,  0^2 >  were 
degraded  based  on  experimental  test  data  for  the  AS/3501-5  system.  Each 
ply  orientation,  for  a  12  in.  panel  radius,  was  evaluated  at  20 
time/ temperature  conditions  that  ranged  from  80  to  300fF,  and  moisture 
concentrations  ranging  from  a  zero  moisture  content  to  Ian  equilibrium 
moisture  distribution.  The  investigation  of  curvature  Iwas  conducted 
only  for  the  (45/-45l2^  laminate  and  at  a  limited  number  of 
time/temperature  conditions.  \ 

The  bifurcation  loads  were  determined  using  the  STAGS-Cl  finite 
elements  shell  analysis  program. c  This  analysis  used  the  prebuckled 
linear  displacement  option  to  calcinate  the  bifurcation  loads.  An 
increase  in  temperature  and  moisture  abso^-pt^on  was  found  to  cause  a 
reduction  in  the  panels  bifurcation  load  rangirtg^rom  a  maximum  of  25.6 
percent  for  the  (0/45/-45/90]g  laminate  to  34.5  percent^for  the 
[45/-45l2s  laminate  for  the  panels  with  a  12  in.  radius ."^This  reduction 
in  the  bifurcation  load  is  significantly  influenced  by  the  change  in 
curvature  at  elevated  temperatures  and  moisture  content. ^The  maximum 
reduction  in  the  bifurcation  load  varied  from  34.5  percenV  for  a  12  in. 
panel  radius  to  11.9  percent  for  a  10,000  in.  panel  radiusXa  decrease  of 
22.6  percent  for  the  l45/-45]2s  laminate.  \ 
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